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Abstract 

Research  was  conducted  to  determine  the  effect  of  F-4E  aircraft 
dynamics  upon  pilot  tracking  capability  in  the  aerial  gunnery  mission. 

An  example  firing  situation  was  postulated  and  the  characteristics  of 
the  resulting  lead-pursuit  trajectory  were  determined.  Mathematical 
representations  of  the  aircraft,  the  flight  control  system  dynamics, 
the  fire  control  system  dynamics,  the  turbulence  environment,  and  the 
pilot  were  used  to  predict  tracking  error.  The  fire  control  system  was 
assumed  to  compute  the  correct  aimpoint,  and  the  projectile  terminal 
miss  distance  due  to  pilot  tracking  error  was  computed  using  an  Iterative 
scheme.  The  results  of  this  analysis  were  then  compared  with  a  baseline 
F-4C  strafing  case  and  found  to  be  similar  in  trend  but  greater  in 
magnitude. 
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ANALYSIS  OF  PILOTED  WEAPON  DELIVERY: 

F-4E  AERIAL  GUNNERY 

I.  Introduction 

The  princ’pal  objective  in  designing  an  aircraft  should  be  the 
creation  of  a  system  which  allows  the  pilot  to  accomplish  the  intended 
mission  with  ma>1mum  effectiveness.  Until  recently,  however,  this  design 
goal  could  not  be  achieved  directly  as  no  analytical  theory  existed  which 
related  aircraft  design  parameters  to  mission  success  criteria. 

The  primary  obstacle  to  developing  such  a  system  theory  of  aircraft 
design  was  the  lack  of  knowledge  of  the  performance  characteristics  of 
the  pilot.  In  many  situations.  Including  the  one  described  In  this 
report,  the  pilot  acts  as  a  feedback  controller  using  information  about 
the  current  aircraft  state  to  generate  commands  to  complete  his  assigned 
mission.  Thus,  the  major  consequence  of  this  inability  to  mathematically 
model  the  pilot  was  that  the  des  gner  could  analytically  define  only  a 
portion  of  the  forward  transfer  function  consisting  of  the  aircraft,  the 
flight  control  system,  the  fire  control  system,  and  any  other  subsystems 
included  In  the  signal  path  between  the  pilot's  input  and  the  aircraft 
controlled  variables.  This  state  of  affairs  produced  two  major  effects. 

First,  aircraft  design  goals  were  stated  in  terms  of  open  loop 
dynamic  parameters  (Ref  4).  Empirical  data  showed  that  if  an  aircraft 
exhibited  certain  response  characteristics,  the  pilot  would  be  able  to 
accomplish  his  mission.  No  claim  was  made  that  the  properly  designed 
aircraft  would  allow  the  pilot  to  do  this  in  the  most  efficient  manner. 

Second,  aircraft  subsystems  were  treated  as  separate  entities 
rather  than  as  components  of  a  total  system.  Since  the  contribution  of 
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each  functional  unit  to  the  dynamic  characteristics  of  the  entire  system 
and  the  Interaction  of  these  units  could  not  be  determined,  the  design 
approach  was  to  ootimlze  each  item  as  a  separate  entity  regardless  of 
the  effect  upon  the  total  system.  Thus,  the  emphasis  of  design  improve¬ 
ment  tended  to  be  concentrated  at  the  subsystem  level  as  evidenced  by 
the  continued  reduction  of  weapon  anomalies,  development  of  more  accurate 
sensors,  improvement  of  computation  schemes,  reduction  of  component  toler¬ 
ances,  and  a  host  of  other  innovations.  It  was  assumed  without  question 
that  these  subsystem  Improvements  would  create  a  more  efficient  system 
(Ref  17:1). 

This  situation  began  to  change  in  the  1960s  when  McRuer,  Levi  son, 
and  many  others,  began  a  series  of  continuing  investigations  into  the 
nature  of  pilot  control  policies  (Refs  9,  10,  and  12).  They  were  able 
to  develop  describing  functions  which  had  the  same  average  output  charac¬ 
teristics  and  produced  the  same  performance  levels  that  a  pilot  did  when 
confronted  with  the  identical  control  task.  These  pilot  describing 
functions  provided  the  designer  with  a  powerful  analytic  tool;  they  could 
be  combined  with  the  aircraft  system  descriptions  to  form  a  closed  loop 
model  of  the  entire  process.  Thus,  the  designer  now  had  the  capability 
to  relate  aircraft  system  parameters  to  mission  performance. 

In  1967,  McRuer  and  Jex  developed  an  analytical  technique  which 
represents  a  systems  approach  to  aircraft  design  (Ref  13:9-28).  Major 
Robert  R.  Rankine,  Jr.,  an  Assistant  Professor  of  Electrical  Engineering 
at  the  Air- Force  Institute  of  Technology,  identified  the  elements  of 
this  theory  for  weapon  delivery  and  analyzed  the  performance  of  the  F-4C 
aircraft  in  the  strafing  role  (Ref  17).  The  basis  for  the  theory  is  the-' 
prediction  of  system  performance  using  an  analytical  model  of  the  pilot; 
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the  method,  which  Is  valid  for  a  wide  variety  of  aircraft  missions, 
proceeds  in  the  following  manner. 

First,  a  n.ssion  is  selected  and  a  performance  criterion  is  formu¬ 
lated  to  measure  the  effectiveness  of  the  pilot-aircraft  system  in 
accomplishing  th- s  mission.  Next,  the  tasks  that  the  pilot  must  perform 
in  order  to  complete  the  mission  are  determined  and  then  all  elements  of 
the  aircraft  and  its  associated  subsystems  which  the  pilot  will  use  to 
accomplish  this  mission  are  identified  and  described  mathematically. 

Models  are  developed  relating  the  aircraft  controlled  variables  to  both 
pilot  and  disturbance  Inputs;  then  mathematical  relationships  are  generated 
which  associate  a  particular  value  of  the  performance  measure  with  each 
aircraft  state  defined  by  the  controlled  variables.  The  end  result  Is 
an  analytical  model  which  relates  the  form  and  values  of  various  parameters 
of  the  pilot-aircraft  system  to  the  measure  of  mission  effectiveness  (Ref  17) 

In  this  study,  the  systems  analysis  method  just  described  is  applied 
to  the  aerial  gunnery  mission.  In  order  to  completely  define  the  situa¬ 
tion,  to  fully  illustrate  the  procedure,  and  to  provide  numerical  results, 
an  example  firing  state  Is  postulated  and  the  F-4E  aircraft  Is  selected 
as  the  system  to  be  analyzed.  To  successfully  complete  this  mission,  the 
pilot  is  required  to  hold  an  aim  point,  the  pipper,  on  the  target  aircraft 
while  firing  the  cannon;  this  task  Is  defined  to  be  target  tracking.  The 
primary  purpose  of  this  thesis  Is  to  show  the  effects  that  flight  control 
system  and  fire  control  system  dynamics  have  upon  pilot  tracking  performance. 

The  analysis  begins  in  Chapter  II  with  a  discussion  of  the  aerial  gun¬ 
nery  mission  and  the  definition  of  a  nominal  firing  state  to  serve  as  an 
illustrative  example  of  this  mission.  Next,  the  pilot  control  tasks  are 
determined  and  the  subsystems  which  are  used  by  the  pilot  to  accomplish 
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the  mission  are  described.  Finally,  an  appropriate  performance  measure 
Is  selected.  In  Chapter  III,  two  mathematical  models  are  developed 
which  relate  pilot  tracking  performance  to  disturbance  Inputs.  This  is 
accomplished  by  perturbing  the  attacking  aircraft  from  Its  nominal 
firing  state  with  turbulence  Inputs  and  then  measuring  the  pilot's 
ability  to  return  the  aircraft  to  the  nominal  state.  A  final  model  Is 
formulated  -In  Chapter  IV  to  assign  a  specific  value  to  the  performance 
measure  for  each  level  of  pilot  tracking  performance.  In  Chapter  V 
the  conclusions  which  can  be  drawn  from  this  analysis  are  presented. 

The  systems  theory  of  aircraft  design  employed  In  this  thesis  1 s  a 
powerful  method  which  has  unique  capabilities  and  great  versatility. 

It  can  be  used  to  predict  the  performance  of  an  aircraft  In  all  weapon 
delivery  roles  and  provides  a  means  for  compering  different  aircraft 
employed  In  the  same  mission.  It  also  allows  the  analytic  evaluation 
of  proposed  design  Improvements.  Since  the  factors  which  degrade  per¬ 
formance  are  analytically  Identified  and  ranked  In  order  of  relative 
Importance  by  this  procedure,  deficiencies  may  be  corrected  In  the  logical 
order  of  their  contribution  to  the  performance  measure.  This  Is  a  powerful 
tool  for  the  manager  as  It  enables  him  to  make  more  effective  use  of  his 
limited  resources.  Applied  on  a  large  scale,  the  method  could  be  used 
to  Identify  basic  airframe  characteristics  and  bese  line  subsystems  In 
the  flight  control  system  and  the  fire  control  system  which  would  be 
required  to  achieve  a  given  acceptable  level  of  mission  performance. 

Finally,  If  probability  of  mission  success  is  used  as  a  performance  cri¬ 
terion,  the  military  strategist  could  employ  this  procedure  to  determine 
the  force  levels  required  to  achieve  given  strategic  and  tactical  objectives 
and  to  compare  his  resources  to  those  of  the  opponent  to  determine  the  proba 
bility  of  success  In  a  given  military  venture  (Refs  18  and  19). 
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1 1 .  The  Aerial  Gunnery  Mission 


Introduction 

In  this  chapter  the  aerial  gunnery  mission  Is  discussed  and  an 
example  firing  situation  is  postulated.  The  closed-loop  control  system 
responsible  for  accomplishing  the  mission  is  described  and  the  functions 
of  various  component's  of  this  system  are  examined.  Finally,  an  appropriate 
performance  criterion  Is  selected. 

Aerial  Gunnery  Mission 

The  objective  of  the  aerial  gunnery  mission  Is  to  score  projectile 
hits  upon  the  target  aircraft;  the  pilot  achieves  this  goal  by  flying 
his  aircraft  to  a  required  attitude  and  holding  this  attitude  while 
firing  the  cannon  (Ref  8:1).  The  fire  control  system  determines  the 
correct  aircraft  attitude  and  supplies  this  information  to  the  pilot  via 
an  aiming  dot.  As  the  attacking  aircraft  moves  In  space,  the  fire  control 
system  continuously  computes  the  firing  states  which  comprise  the  lead- 
pursuit  trajectory.  The  pilot  remains  on  this  trajectory  by  superimposing 
the  aiming  dot  on  the  target  aircraft;  a  projectile  hit  Is  guaranteed  if 
the  pilot  fires  the  cannon  while  on  this  nominal  trajectory.  There  are 
literally  an  infinite  number  of  lead-pursuit  trajectories  based  upon  the 
target  and  the  attacking  aircraft  states;  therefore,  an  example  firing 
situation  will  be  stated  and  used  to  highlight  certain  aspects  of  the 
aerial  gunnery  mission  and  to  Illustrate  the  analytical  techniques  applied 
In  this  study. 

Example  Firing  State  __ 

The  description  of  the  example  firing  state  begins  with  the  definition 
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of  two  terms.  The  weapon  1 Ine  Is  defined  to  be  the  extension  of  the 
cannon  bore  and  !s  col  Inear  with  the  muzzle  velocity  vector  of  the  pro¬ 
jectile;  the  tra:king  line  Is  the  line  from  the  pilot's  eye  which  passes 
through  the  aiming  dot.  The  weapon  line  and  the  tracking  line  are  ini¬ 
tially  considered  to  be  col  inear;  then,  the  target  aircraft  is  located 
at  a  range  of  2000  feet  along  the  combined  weapon  and  tracking  lines. 

In  the  absence  of  aircraft  motion  and  field  effects  (gravity,  projectile 
drag),  firing  the  aircraft  cannon  from  this  position  would  result  In  a 
hit  on  the  target. 

Next,  the  two  aircraft  are  set  in  motion.  Initially,  the  attacking 
aircraft's  velocity  Is  0.9  Mach  and  is  assumed  to  be  along  the  combined 
weapon  and  tracking  lines.  For  computational  purposes.  It  is  necessary 
to  reference  the  target  aircraft's  velocity  vector  to  the  tracking  line. 
Hence,  a  tracking  line  reference  system  is  defined  and  illustrated  in 
Fig.  1,  where 

TLX  is  measured  along  the  tracking  line, 

TLY  Is  measured  In  the  plane  perpendicular  to  the  tracking 
line  and  parallel  to  the  local  horizontal, 

TLZ  Is  measured  In  the  plane  .perpendicular  to  the  tracking 
line  and  Is  also  perpendicular  to  the  TLY  axis, 
p  Is  the  target  aircraft  heading  measured  positively  In 
a  clockwise  sense  about  the  TLZ  axis,  and 
X  Is  the  target  aircraft  flight  path  angle  measured  posi¬ 
tively  in  a  counterclockwise  sense  about  the  displaced 
TLY  axis. 

In  this  example,  the  target  aircraft  has  a  velocity  of  .8  Mach,  a  heading^ 
of  0°,  and  a  flight  path  angle  of  0°;  it  is  further  assumed  that  the  target 
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Is  in  tmaccelcrated  flight  and  will  remain  in  that  condition  during  the 
projectile  time-of-fl ight. 

When  field  effects  and  aircraft  motion  are  included  in  the  example 
firing  situation.  It  becomes  necessary  to  create  the  proper  angle  between 
the  weapon  line  and  the  tracking  line  while  holding  the  tracking  line  on 
the  target  aircraft  to  Insure  projectile  hits.  This  angle  is  called  the 
prediction  angle  and  is  determined  by  the  fire  control  system  of  the 
attacking  aircraft.  The  calculation  of  the  prediction  angle  for  this 
example  firing  situation  is  covered  In  detail  in  Chapter  IV;  additionally, 
some  of  the  functional  aspects  of  the  computation  of  this  angle  by  the 
fire  control  system  are  discussed  in  the  next  section. 

To  complete  the  specification  of  the  example  firing  state,  both  air¬ 
craft  are  assumed  to  be  flying  at  an  altitude  of  15,000  feet  mean  sea 
level  and  the  attacking  aircraft,  an  F-4E,  is  operating  at  a  weight  of 
46,600  pounds. 
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This  case  represents  one  of  the  simplest  firing  situations  that  could 
occur  for  the  attacking  aircraft.  Both  aircraft  are  travelling  in  the 
same  direction,  the  F-4E  is  within  firing  range,  and  is  overtaking  the 
target  aircraft  directly  from  the  rear.  The  parameters  of  the  firing 
state  were  dictated  by  the  limits  of  validity  of  the  pilot  model  and  by 
the  desire  to  reduce  computational -errors  within  the  fire  control  system 
to  zero.  In  spite  of  the  seeming  simplicity  cf  the  situation,  however, 
the  aircraft  dynamics  do  exert  a  powerful  influence  upon  the  pilot  track¬ 
ing  performance  as  will  be  shown  and,  since  the  objective  of  this  work 
is  to  demonstrate  the  relationship  between  these  two  items,  the  firing 
situation  is  sufficiently  rigorous  for  the  purposes  of  this  analysis. 

A  summary  of  the  parameters  of  the  example  firing  state  is  presented 
in  Table  I. 

TABLE  I 

Firing  State  Parameters 


Target  range 

2000  ft 

Target  aircraft: 

Heading 

0  deg 

Flight  path  angle 

0  deg 

Mach  number 

0.8 

Velocity 

845.9  ft/sec 

Attacking  aircraft: 

Mach  number 

0.9 

Velocity 

951.6  ft/sec 

Altitude 

15,000  ft 

Weight 

46,159  lb 
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Aerial  Gunnery  System 

General  Discussion.  The  aerial  gunnery  system  Is  defined  to  en¬ 
compass  the  pllct,  the  aircraft  and  its  associated  dynamics,  the  flight 
control  system  including  all  stability  augmentation,  the  fire  control 
system  including  all  appropriate  measurements,  computations,  and  displays, 
and  the  aircraft  mounted  cannon,  a  20  mm  Gatling  gun. 

For  purposes  of  this  analysis,  it  Is  assumed  that  the  attacking 
aircraft  has  achieved  stabilized  tracking  on  the  nominal  lead-pursuit 
trajectory.  Thus,  the  radar  Is  locked  on  to  the  target,  the  fire  control 
system  Is  computing  and  displaying  the  correct  aimpcint,  and  the  pilot 
has  placed  the  aimpoint  on  the  target  aircraft.  The  net  effect  Is  that 


there  are  no  transients  within  the  aerial  gunnery  system. 

Several  salient  features  of  the  aerial  gunnery  system  will  now  be 

* 

discussed;  the  Intent  is  to  provide  insight  into  the  operation  of  the 
system,  derive  the  control  laws  which  must  be  Implemented  if  the  mission 


is  to  succeed,  and  to  describe  the  limitations  of  the  system. 


Fire  Control  System.  The  F-4E  Fire  Control  System  consists  of  the 
sight  reticle,  the  lead-computing  amplifier,  and  the  lead-computing 


gyroscope;  the  purpose  of  these  components  Is  to  solve  Eq  (1)  for  the 
prediction  angle  and  display  this  Information  to  the  pilot'  (Ref  7:2-1). 
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where 

X  =  vector  prediction  angle, 

Tf  =  projectile  time  of  flight, 

I  *  angular  velocity  of  the  llne-of-slght, 

Jv  *  ballistic  jump  parameter, 

Va  ■  attacker  velocity, 

Vf  ^./vricasc. project i; la. velocity  relative  to  the  attacking  aircraft, 
a  =  gun  angle  of  attack, 

5  c  unit  vector  along  tracking  line, 
an  ■  normal  acceleration, 

D  *  range, 

Ag  *  drag  coupling  coefficient,  and 
z  *  unit  vector  along  vertical. 

Unfortunately,  the  system  can  compute  only  the  first  three  terms  of  Eq  (1) 
due  to  measurement  problems  (Ref  7:2-2);  Therefore,  If  there  are  to  be 
no  computational  errors  In  the  probl'em,  the  example  firing  situation  must 
be  constructed  in  such  a  manner  that  the  final  two  terms  are  truly  zero; 
it  Is  for  this  reason  that  both  the  target  aircraft  and  the  attacking 
aircraft  were  constrained  to  unaccelerated  flight.  Unfortunately,  this 
requirement  also  reduces  the  first  term,  lead  for  target  motion,  to  zero; 
hence,  the  final  computations  include  only  gravity  drop,  jump  angle,  and 
projectile  drag  effects.  For  the  example  firing  situation,  then,  Eq  (1) 
reduces  to  Eq  (123),  Eq  (128),  and  Eq  (130)  of  Chapter  IV  and  the  pre¬ 
diction  angle  calculated  by  the  fire  control  system  is  the  same  as  the 
true  prediction  angle. 

The  lead  pursuit  tracking  task  is  an  Inherently  unstable  situation^ 
to  remedy  this  fact,  a  lag  Is  incorporated  in  the  fire  control  system. 
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This  lag  Is  called  the  "stability  number"  and  Is  .25  for  the  F-4E  (Ref 
8:15).  Thus,  there  Is  a  delay  between  the  calculation  of  a  given  pre¬ 
diction  angle  and  the  display  of  this  Information  to  the  pilot.  In 
certain  dynamic  situations,  then,  the  sight  reticle  Is  Indicatl  ig  an 
incorrect  aiming  point.  Since  sight  instability  Is  most  pronounced  during 
the  initiation  of  target  tracking  and  under  high  attacking  aircraft  accel¬ 
eration  conditions  and,  since  the  example  firing  situation  is  an  extremely 
benign  case,  the  sight  stability  number  will  be  considered  to  be  zero  In 
this  study  and  It  will  be  assumed  that  the  aim  dot  Is  indicating  the 
correct  firing  information  at  every  instant  of  time. 

This  completes  the  discussion  of  the  unique  features  and  limitations 
of  the  fire  control  system.  There  are  several  other  error  generating 
factors  In  this  situation;  they  all  have  a  relatively  minor  effect  on  the 
problem  solution,  are  briefly  discussed  in  Chapter  IV,  and  are  neglected 
for  the  purposes  of  this  analysis.  Fire  control  system  computational 
errors  are  ignored  for  the  remainder  of  this  study;  if  the  pilot  fires 
the  cannon  with  the  aim  dot  on  the  target,  it  is  assumed  that  the  pro¬ 
jectile  will  hit  the  target  aircraft. 

Sight  Reticle.  The  information  necessary  to  create  the  required 
prediction  angle  is  presented  to  the  pilot  in  the  form  of  a  sight  reticle; 
this  reticle  is  a  visual  pattern  displayed,  on  the  optical  sight  unit. 

The  optical  sight  unit  is  mounted  on  top  of  the  center  of  the  instrument 
panel  and  is  the  only  instrument  used  by  the  pilot  to  accomplish  the 
aerial  gunnery  mission.  The  sight  reticle  is  illustrated  in  Fig.  2 
(Ref  20:1-121). 

The  sight  reticle  is  initially  aligned  with  the  weapon  line  as 
discussed  in  the  example  firing  situation.  After  a  radar  lock-on  has  been 
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achieved  and  the  prediction  angle  has  been  computed,  the  sight  reticle 
moves  downward  on  the  optical  sight  unit  through  the  prediction  angle; 
the  pilot  then  flies  the  pipper  (aim  dot)  back  up  to  the  target  aircraft. 
This  action  raises  the  weapon  line  above  the  target  in  the  target's  direc¬ 
tion  of  travel  by  the  required  prediction  angle;  thus,  the  firing  attitude 
required  to  achieve  a  hit  on  the  target  has  beei  established. 

The  fire  control  system  will  generate  prediction  angles  for  ranges 
varying  from  900  feet  to  4000  feet.  The  pilot  may  fire  then  anytime  the 
range  bar  indicates  a  range  between  these  two  extremes  (Ref  20:124).  The 
range  bar  indicates  a  range  of  6700  feet  when  fully  expanded,  4000  feet 
at  the  state  shown  In  Fig.  2,  and  900  feet  when  fully  contracted  at  the 
bottom  of  the  reticle.  The  bar  operates  like  a  thermometer,  reducing 
the  amount  of  red  field  displayed  as  range  decreases.  The  gun  Is  opti¬ 
mized  for  a  range  of  2250  feet  and,  without  radar  range  information,  the 
prediction  angle  Is  calculated  on  the  basis  of  a  1500-foot  range.  The 
pilot  Is  therefore  more  likely  to  fire  In  this  interval  of  ranges  and,  as 
a  result,  a  target  range  of  2000  feet  is  used  in  the  example  firing  situa¬ 
tion.  There  Is  one  other  set  of  competing  considerations  In  selecting  a 
firing  range.  Changes  in  the  target  acceleration  during  the  projectile 
time  of  flight  will  result  in  a  terminal  error.  This  error  can  be  mini¬ 
mized  by  firing  at  shorter  ranges;  however,  the  attacking  aircraft  runs 
the  risk  of  colliding  with  the  target  or  hitting  target  debris.  All  of 
these  factors,  then,  generally  result  in  firing  ranges  from  1500  feet  to 
2500  feet  and,  therefore,  the  firing  range  in  the  example  firing  situation 
Is  certainly  realistic. 

Pilot  Control  Tasks.  As  the  attacking  aircraft  moves  along  the  nominal 
lead-pursuit  trajectory,  it  is  disturbed  from  this  path  by  atmospheric 
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turbulence;  the  pilot  views  this  perturbation  from  the  nominal  state  as 
both  lateral  and  longitudinal  displacement  of  the  plpper  from  the  target 
aircraft.  !r  order  to  accomplish  the  aerial  gunnery  mission,  the  pilot 
must  return  the  plpper  to  the  target;  to  do  this,  he  must  perfonn  three 
distinct  control  tasks. 

The  pilot's  first  concern  Is  to  maintain  the  nominal  trajectory  to 
insure  .future  firing  opportunities  and  he  achieves  this  by  properly  con¬ 
trolling  the  roll  angle  of  the  attacking  aircraft.  Referring  to  Fig.  1, 
page  7,  the  TLY-TLZ  plane  Is  called  the  maneuver  plane.  The  attacking 
aircraft  must  move  In  the  same  direction  in  this  plane  as  the  target 
aircraft  does  to  insure  continued  tracking.  The  pilot  accomplishes  this 
goal  by  predicting  the  future  position  of  the  target  from  Its  present 
attitude,  velocity  and  acceleration  and  then  flying  the  plpper  to  this 
predicted  position.  This  estimate  is  easy  to  make  in  the  example  firing 
situation;  the  target  aircraft  is  located  directly  ahead  of  the  attacker, 
has  no  velocity  in  the  maneuver  plane,  and,  as  It  can  only  accelerate  in 
the  maneuver  plane  along  its  lift  line,  its  vertical  tail  points  in  the 
direction  of  its  next  most  probable  position.  Therefore,  the  pilot  will 
attempt  to  maintain  the  same  roll  attitude  as  the  target  aircraft;  this 
will  enable  him  to  maintain  a  stable  lead-pursuit  trajectory  by  using 
his  acceleration  capability  to  move  with  the  target  in  the  maneuver  plane. 

The  example  firing  situation  contains  one  subtle  element  which  is 
crucial  to  the  success  of  this  analysis.  For  the  analytic  pilot  modelling 
technique  to  be  valid,  the  pilot  must  have  an  explicit  error  signal  upon 
which  to  operate  (Ref  12:240).  Here,  the  upper  index  of  the  sight  reticle 
and  the  vertical  tail  of  the  target  aircraft  provide  an  explicit  measure 
of  bank  angle  error.  As  the  target  aircraft  achieves  an  angle  off,  however. 
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It  begins  to  have  a  velocity  component  In  the  maneuver  plane,  and  bank 
angle  no  longer  provides  an  adequate  cue  to  the  future  position  of  the 
target.  The  pilot  then  relies  on  his  experience  to  determine  the  probable 
future  position  of  the  target  from  the  complex  mixture  of  target  attitude, 
velocity,  and  acceleration.  Since  the  pilot  Is  no  longer  operating  on  an 
explicit  error  signal,  one  of  the  constraints  of  the  analytical  pilot 
model  ling  technlcje  is  no  longer  satisfied  and  there  is  no  assurance  that 
the  model  can  be  ised.  Thus,  the  parameters  of  the  example  firing  situa¬ 
tion  were  selected  to  specifically  satisfy  the  requirements  of  the 
analytical  pilot  model. 

Having  achieved  a  stable  nominal  trajectory,  the  pilot  next  acts  to 
eliminate  the  lateral  displacement  of  the  pipper.  The  lateral  error  in 
plpper  placement  perceived  by  the  pilot  is  due  to  yaw  perturbation  M  and 
a  component  of  roll  ($)  caused  by  the  depression  of  the  tracking  line 
below  the  velocity  vector,  the  axis  of  roll  of  the  aircraft;  the  geometry 
of  the  situation  Is  shown  in  stability  axis  in  Fig.  3.  From  Fig.  3  and 
the  relations  for  spherical  triangles, 

sin  <p  _  sin(Ti/2)  _  sin(7r/2  - $)  ^2) 

sin  C  sin(I-a0)  sin  D 

where  I  is  the  sight  depression  angle  in  radians,  and  ag  is  the  aircraft 
angle  of  attack  on  the  nominal  trajectory.  Since  sin(r/2)=l  and 
sin(Ti/2  - <p)  =  cos  $,  this  gives  the  relations: 

sin  C  =  sin  $  sin(Z-aQ)  (3) 

sin  D  =  cos  $  sin(I  -aQ)  (4) 

When  I  >  <50 »  a  positive  roll  perturbation  causes  the  pipper  to  mcve  to  the 
left  of  the  target  while  a  positive  yaw  perturbation  always  causes  the 
pipper  to  move  to  the  right  of  the  target.  Assuming  the  longitudinal 
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ig.  3.  The  Geometry  of  Lateral  Pipper  Displacement 
(Adapted  from  Ref  17:25) 
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variables  to  be  at  their  nominal  values,  the  lateral  tracking  error  ob¬ 
served  by  the  p*.  1  ot  can  be  expressed  as 

vi  =  tfi  -  C  (5) 

where  from  Eq  (3), 

C  =  sin-1  [sin  $  sin(Z  -o0)]  (6) 

Consistent  with  the  small  perturbation  analysis  being  accomplished  in  this 
study,  Eq  (5)  ma>  be  simplified  to 

v  =  -  (Z  -  a0>4>  (7) 

Therefore,  the  lateral  tracking  task  consists  of  controlling  a  linear 
combination  of  perturbation  yaw  and  roll  (Ref  17:22-24). 

Finally,  the  pilot  must  eliminate  the  longitudinal  displacement  of 
the  pipper  from  the  target  aircraft.  When  the  attacking  aircraft  is 
disturbed  from  the  nominal  trajectory,  a  pitch  perturbation  angle  (e) 

Is  created.  The  lead  computing  gyroscope  senses  the  rate  of  this  angle 
and  commands  a  new  pipper  position  on  the  optical  sight  unit.  Thus,  the 
pilot  observes  and  must  control  a  combination  of  perturbation  pitch  angle 
(e)  and  perturbation  pipper  position  (ep). 

Control  Input.  The  final  function  in  the  sequential  description  of 
the  aerial  gunnery  system  Is  the  generation  of  the  appropriate  controls 
to  place  the  pipper  on  the  target.  In  the  longitudinal  tracking  case, 
a  pitch  input  at  the  stick  is  the  only  control  that  is  available  to  the 
pilot.  In  the  roll  angle  and  lateral  tracking  tasks,  however,  the  pilot 
has  two  control  inputs  available;  the  pilot  may  command  aileron,  rudder, 
or  a  combination  of  these  two  controls. 

The  pilot  can  adopt  three  distinctly  different  strategies  in  the  case"* 
of  lateral-directional  control  inputs.  He  may  use  aileron  or  rudder  alone 
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or  he  may  use  one  Input  as  a  primary  control  and  the  other  Input  as  a 
secondary  control.  An  example  of  this  policy  Is  the  Introduction  of 
aileron  to  create  a  given  bank  angle  and  turn  rate  and  then,  In  a  secon¬ 
dary  role,  the  introduction  of  rudder  to  prevent  slipping  or  skidding. 

Here  the  pilot  uses  aileron  to  accomplish  his  task  and  rudder  to  reduce 
undesirable  consequences  of  the  primary  control  action.  Finally,  a  pilot 
may  use  both  rudder  and  aileron  in  a  primary  control  sense.  An  excellent 
example  of  this  s* tuatlon  is  the  simultaneous  application  of  rudder  and 
aileron  to  create  a  sideslip  for  a  crosswind  landing.  Once  the  controls 
have  taken  effect,  the  pilot  will  alter  one  or  both  inputs  to  achieve  his 
objective.  Here,  one  is  unable  to  describe  either  input  as  a  secondary 
control.  A  well -trained  pilot  will  select  the  optimum  control  or  mix  of 
controls  to  accomplish  his  mission;  further,  the  pilot  with  experience 
In  a  particular  aircraft  has  learned  exactly  the  magnitude  and  phasing  of 
control  Inputs  required  to  successfully  complete  his  mission. 

In  this  analysis,  an  aileron  command  Is  the  control  Input  applied  to 
accomplish  the  roll  angle  and  lateral  tracking  tasks.  The  author's  exper¬ 
ience  in  tactical  aircraft  similar  to  the  F-4E  supports  this  assumption 
for  the  low  angle  of  attack  situation  studied  here.  Since  the  author  has 
no  flying  experience  In  the  F-4E  and  therefore  is  unfamiliar  with  the  pilot 
control  strategies  peculiar  to  that  aircraft,  no  secondary  rudder  commands 
will  be  used  in  the  example  firing  situation  to  counter  the  effects  of  pro- 
verse  yaw.  However,  at  high  angles  of  attack  and  at  high  normal  acceleration 
levels,  the  rudder  becomes  an  extremely  effective  control  for  maneuvering 
and  tracking.  Therefore,  in  each  firing  situation,  some  thought  must  be 
given  to  pilot  control  policies;  It  cannot  be  automatically  assumed  that 
the  pilot  will  only  use  aileron  commands  to  accomplish  the  aerial  gunnery 
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missions.  The  Intent  of  this  discussion,  then,  Is  to  present  an  Insight 
Into  the  variety  of  control  options  open  to  the  pilot  In  order  that  the 
reader  may  more  fully  understand  the  complexities  of  this  particular 
analysis. 

Performance  Criterion 

A  performance  criterion  must  be  -selected  to  measure  the  pilot's 
effectiveness  In  accomplishing  the  aerial  gunnery  mission.  The  magni¬ 
tude  of  the  perturbations  of  the  controlled  variables,  the  projectile 
terminal  miss  distance,  and  probability  of  kill  are  all  meaningful 
expressions  of  the  degree  of  task  completion. 

In  this  study,  the  magnitude  of  the  perturbations  of  the  controlled 
variables  Is  selected  as  the  performance  criterion;  the  form  in  which 
these  magnitudes  are  expressed,  root  mean  square  intensity,  Is  explained 
in  Chapter  III.  This  criterion  was  selected  as  it  provides  a  non-dimensional 
n.jasure  which  can  easily  be  compared  to  similar  quantities  in  other  flight 
conditions  or  missions.  Miss  distance,  on  the  other  hand,  is  meaningful 
only  If  target  range  is  constant;  one  key  variable  in  the  aerial  gunnery 
task  is  the  target  range  and,  consequently,  miss  distance  constrains  the 
analysis  if  It  Is  used  as  a  performance  measure.  Probability  of  kill 
places  emphasis  upon  probability  of  mission  accomplishment  and  therefore 
Is  very  useful  to  the  planner.  However,  it  does  not  exhibit  the  source 
or  magnitude  of  the  Individual  errors  and,  thus,  is  not  a  useful  tool 
for  the  engineer. 
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III.  Prediction  of  Tracking  Performance 

The  purpose  of  this  chapter  Is  to  describe  In  detail  the  develop¬ 
ment  of  two  mathenatical  models  relating  gust  Inputs  to  changes  in  the 
controlled  variables  of  the  attacking  aircraft.  The  great  significance 
of  these  equations  Is  that  they  will  permit  the  quantitative  evaluation 
of  perturbations  In  the  tracking  variables  caused  by  a  given  turbulence 
input.  Thus,  the  nodels  will  allow  prediction  of  pilot  tracking  performance. 

The  critical  element  In  the  formulation  of  these  relationships  Is  the 
representation  of  the  pilot  as  a  linear  element  In  a  complex  control  sys¬ 
tem.  Consequently,  prior  to  deriving  the  desired  transfe  unctions,  the 
mathematical  modeling  of  the  pilot  will  be  discussed. 

Analytic  Pilot  Modeling 

The  pilot  Is  a  complex,  non-linear,  time  varying,  controller  who 
defies  complete  mathematical  description;  however.  In  certain  limited,  well 
defined  situations,  a  pilot's  control  policy  is  so  consistent  that  It  can 
be  represented  by  a  linear  mathematical  model.  The  model  Is  not  constrained 
to  correspond  In  structure  or  methodology  to  the  pilot  but  must  exhibit, 
on  the  average,  the  same  output  characteristics  for  a  given  input  and 
achieve  the  same  measure  of  performance  as  a  pilot  would  when  confronted 
with  the  Identical  task. 

One  situation  for  which  a  well -documented  and  widely  used  pilot 
model  exists  Is  the  compensatory  tracking  task.  In  this  situation,  the 
pilot  observes  the  difference  between  a  stationary  random  input  and  the 
output  of  a  linear  time-invariant  process  and  uses  the  perceived  error  to  ^ 
generate  a  manual  control  Input  to  the  process  (Ref  17:68-71). 
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Examination  of  the  firing  situation  postulated  In  the  previous 
chapter  reveals  that  It  exhibits  all  the  characteristics  of  the  compen¬ 
satory  tracking  task.  First,  the  pilot  Is  controlling  an  aircraft  and 
Its  associated  dynamics  which  can  be  accurately  represented  as  a  linear, 
time-invariant  s>stem.  Next,  this  plant  Is  disturbed  from  Its  nominal 
trajectory  by  a  stationary  random  Input  In  the  form  of  atmospheric  tur¬ 
bulence.  The  ptlot  views  this  perturbation  from  the  nominal  state  as  a 
displacement  of  the  pipper  from  the  target  aircraft  and  he  generates  a 
control  Input  In  the  form  of  a  stick  or  rudder  movement  to  again  super¬ 
impose  the  pipper  on  the  target.  Thus  all  the  requirements  of  the 
compensatory  tracking  task  have  been  satisfied  by  the  example  firing  situa¬ 
tion  and  It  Is  valid,  for  purposes  of  analysis,  to  represent  the  pilot  as 
a  linear  element  In  the  control  system. 

The  model  of  the  compensatory  tracking  task  Is  Illustrated  In 
Fig.  4.  The  pilot  Is  represented  by  a  linear  element  and  a  remnant,  the 
latter  accounting  for  that  part  of  the  pilot's  output  which  Is  not  linearly 
correlated  with  the  Input.  The  linear  element  of  the  model  must  account 
for  the  adaption  of  the  pilot  and  consists  of  two  distinct  parts:  (1)  a 
generalized  describing  function  form  representing  the  pilot's  capability 
to  compensate,  and  (2)  a  set  of  rules  for  adjustment  of  the  parameters  of 
the  describing  function  for  the  Input  and  plant  characteristics.  The 
model  permits  prediction  of  pilot  performance  In  terms  of  root  mean  square 
error  (oe)  when  the  bandwidth  and  root  mean  square  amplitude  of  the  Input 
(o$)  and  the  aircraft  dynamics  are  known  (Ref  17:68-69). 

A  simple  describing  function  form  that  can  be  used  to  represent  the 
pilot  in  the  compensatory  tracking  task  Is  (Ref  12:234) 

1  +  ju)TLN 


3-]WXe 


(8) 
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Fig.  4.  Quasl-Llnear  Model  of  Compensatory  Tracking  by  a  Pilot 
(Adapted  from  Ref  12:240) 


where 

Kp  =  pilot  static  gain, 

Tl  =  lead  time  constant, 

Tj  =  lag  time  constant,  and 

xe  =  effective  time  delay  Including  transport  delays  and  high 
frequency  neuromuscular  lags. 

The  rules  for  adjusting  the  parameters  of  Eq  (  8)  may  be  summarized  as 
an  attempt  to  achieve  an  open-loop  transfer  function  near  the  crossover 
frequency  of  the  form 


yy  =!£e-jurre  (9) 

pc  jw 

where 

Yp  =  pilot  describing  function  from  Eq  (8  ), 

Yc  =  transfer  function  of  the  controlled  element, 
uc  =  crossover  frequency  (loop  gain),  and 
te  =  effective  time  delay. 

O 


22 


GGC/EE/72-3 


The  transfer  function  of  the  controlled  element  Yc  must  Include  the 
dynamics  of. the  aircraft,  the  flight  control  system,  the  fire  control 
system,  and  the  appropriate  displays.  A  comprehensive  summary  of  the 
rules  which  establish  the  parameters  of  Eq  (8)  for  both  single  axis 
and  multi-axis  compensatory  tracking  can  be  fourd  in  the  paper  by  McRuer 
(Ref  12:231-248). 

In  the  example  firing  situation,  the  pilot  must  perform  three  com¬ 
pensatory  tracking  tasks.  In  each  case,  an  analytic  pilot  model  of  the 
form  given  in  Eq  (8)  is  employed  and  the  adaption  rules  are  discussed 
in  greater  detail  and  applied  to  determine  the  parameters  of  the  model  for 
that  particular  task. 

Lateral  Input  and  Controlled  Element 

In  the- example  firing  situation,  the  pilot  has  established  himself 
on  a  nominal  lead-pursuit  trajectory  and  is  presented  with  a  series  of 
firing  opportunities  which,  if  executed,  would  result  in  projectile  hits 
upon  the  target  aircraft.  However,  as  he  moves  along  this  nominal  trajec¬ 
tory,  his  aircraft  is  subject  to  two  lateral  disturbance  inputs,  a  vertical 
turbulence  spectrum  which  induces  a  roll  rate  and  a  horizontal  turbulence 
spectrum  which  produces  both  a  sideslip  and  a  yaw  rate.  These  gust  inputs 
cause  the  aircraft  to  be  perturbed  off  the  nominal  trajectory  and  the 
pilot  observes  this  perturbation  as  a  lateral  displacement  of  the  pipper 
from  the  target. 

The  roll  rate  disturbance  produced  by  the  continuous  random  vertical 
turbulence  spectrum  can  be  mathematically  represented  by  a  low-pass  filtered 
white  noise  source.  The  filter  dynamics  can  be  represented  as  shown  in 
Fig.  5,  where 
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Fig.  5.  A  Moeel  of  the  Effect  of  Vertical  Susts  on  the  Lateral 
Dynamics  of  an  Aircraft 


Tpg(s)  =  V(1/WJ0) 


VC.  8(^/4^  X/3]\ 

- - - )aw 

1  +  4bs/irUo  / 


n2  =  white  noise  signal,  and 

crw  -  rms  intensity  of  the  vertical  gust  in  ft/sec 


for  the  example  firing  situation: 

=  1750  ft 
b  =  38.41  ft 


=  scale  factor, 

=  aircraft  wing  span,  and 


Hence, 


Uq  =  951.63  ft/sec  =  aircraft  speed. 


P£_  C-0245)ow 


nz  (s  + 19.46) 

Similarly,  the  sideslip  and  yaw  rate  disturbances  created  by  the 
continuous  random  horizontal  turbulence  spectrum  can  be  represented  by  a 
low-pass-filtered  white  noise  source  which  is  uncorrelated  with  the  noise 
source  producing  the  rolling  perturbation.  The  filter  dynamics  can  be 
modelled  as  shown  in  Fig.  6,  where  (Ref  4:459) 


Tv  (s)  =  av^Ly/nUo) 


/l  +  V(3)Lvs/U0' 
l  (1  +  Lys/Uo)2 


m  =  a  white  noise  signal, 
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Fig.  6.  A  Model  of  the  Effect  of  Side  Gusts  on  the  Lateral 
Dynamics  of  an  Aircraft 


o  =  rms  intensity  of  the  side  gust  in  ft/sec. 


\  =  145b1 /3 

and 

_  R<s>  -  (1/U°)s 
uo  1  +  (3b/wU0)s 

for  the  example  firing  situation: 


Hence, 


h  =  15,000  ft  =  aircraft  altitude,  and 
Ly  =  3576.01  ft  =  a  scale  factor. 


6  -.000530(3 +  .15)o 

_£  = _ 1 

ni  (S  +  .27)2 


r  -(.0137)(s2  +  .15s)o 
= _ v 

ni  (s + 25.94)(s  +  .27)2 


The  lateral -directional  equations  of  motion  of  the  aircraft  represent 
a  portion  of  the  controlled  element  dynamics  as  well  as  additional  filtering 
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for  the  two  gust  Inputs.  For  the  gust  inputs,  as  well  as  rudder  and 
aileron  Inputs,  the  lateral  directional  equations  are  (Ref  1:112) 
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(17) 


The  aircraft  dynamics  are  altered  by  the  stability  augmentation 

systems  incorporated  in  the  roll  and  yaw  control  channels.  Hence  they 

must  be  included  in  the  above  set  of  matrix  equations. 

The  roll-axis  control  system  can  be  depicted  as  (Ref  17:156-157) 

<5a(s)  -  A(s)[6pA(s)  -  B(s)p(s)]  (18) 

0 

where 

6  (s)  =  aileron  deflection, 

u 

A(s)  =  power  cylinder  dynamics, 

6pA(s)  =  lateral  control  stick  deflection,  and 
B(s)  =  roll  rate  feedback  dynamics. 

Similarly,  the  yaw  axis  control  system  can  be  represented  as 

6r(s)  =  C(s)[5pp(s)  +  D(s)r(s)  +  E(o)3(s)]  -  (19) 

where 

«r(s)  =  rudder  deflection  (positive  left), 

C(s)  =  power  cylinder  dynamics  and  rudder  flexibility, 

6pR(s)  =  rudder  pedal  deflection  (positive  left  forward), 

D(s)  =  yaw-rate  feedback  dynamics,  and 
E(s)  =  sideslip  feedback  dynamics. 
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Substituting  Eqs  (18)  and  (19)  Into  Eq  (17)  gives  the  following  matrix 
equation  relating  the  aircraft  dynamic  response  to  pilot  and  turbulence 
Inputs: 


s-Lp+ABL5a 

-Np+ABN6a 

-sYp-Y*+ABY6fls 


-Lr-CDL5r  -Lg->CEL5r  p 

s-Nr-CDN6r  -NB-CEN5r  r 

-sYp-Y^-CDY^s  s2-sYrCEY5l,s  B 


s-Np-CDN* 
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Np 

CY*rs 
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sY 

r 

Le  6prI 
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•b  n 


The  above  matrix  equation  is  a  lateral -directional  mathematical  model 
of  the  attacking  aircraft  which  the  pilot  must  control  to  successfully 
track  the  target.  Cramer's  rule  (Ref  23:32)  can  be  applied  to  Eq  (20) 
to  obtain  the  transfer  functions  relating  the  aircraft  controlled  vari¬ 
ables  and  s  to  both  gust  and  pilot  inputs.  At  this  point,  then, 
these  equations  can  be  used  to  quantitively  evaluate  the  perturbations  in 
<fr,  <j>,  and  3  caused  by  the  turbulence  input  and  also  to  determine  the  pilot 
actions  necessary  to  restore  the  aircraft  to  the  nominal  lead-pursuit  tra¬ 
jectory.  The  numerical  entries  for  each  element  of  this  matrix  equation 
are  listed  in  Appendix  A,  and  will  be  substituted  in  the  appropriate  posi¬ 
tions  in  the  succeeding  equations. 

Roll  Control  ^ 

At  this  point,  the  attacking  aircraft  has  been  perturbed  off  its 
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nominal  lead-pursuit  trajectory  and  the  pilot  must  generate  an  appropriate 
control  to  bring  the  pipper  back  on  the  target  and  consequently  to  return 
the  aircraft  to  the  nominal  flight  path.  As  discussed  in  the  preceding 
chapter,  the  pilot's  first  act  is  to  reduce  roll  angle  error  to  zero  by 
compensatory  tracking  in  order  to  insure  a  stable  lead-pursuit  trajec¬ 
tory.  Therefore,  an  analytic  pilot  model  must  be  developed,  included  in 
the  system  modei ,  and  a  loop  closure  must  be  made  to  describe  the  pilot's 
control  of  roll  angle  perturbations. 

The  effective  time  delay  parameter  of  Eq  (0)  is  determined  by  the 
characteristics  of  the  disturbance  input.  Using  Cramer's  rule,  the  fol¬ 
lowing  transfer  functions  are  derived: 

.024(s) (s  +  .15)(s  +  .50)(s  +  3. 84)(s- 3.88)o 

■i-i - v- -  (21) 

ni  (s  +  .19)(s  +  .27)2(s2  +  1.77s  +  1.14)(s2  +  11.67s  +  52.80) 

-.0242(s  +  ,80)(s  +1.21)o 

A.  =  - - -  (22) 

nz  (s  +  ,19)(s2  +  1.77s +  1.14)(s2  +  11.67s +  52.80) 

Further  (Ref  4:435), 

0^  =  V(Ly/L^)o^  =  1.43o^7  (23) 

where  is  given  by  Eq  (13),  and  is  given  following  Eq  (10). 

The  symbol,=,  indicates  that  terms  with  a  frequency  above  20  radians/ 
second  have  been  replaced  in  the  transfer  function  by  a  gain.  Computations 
showed  that  this  procedure  altered  the  root  mean  square  tracking  variable 
errors,  the  primary  measure  of  performance  in  this  analysis,  in  the  second 
or  third  significant  digit  of  computation.  Hence,  little  accuracy  was 
sacrificed  but  the  computations  were  simplified  considerably  by  this 
procedure. 

A  log-magnitude  plot.  Fig.  7,  shows  that  the  dominant  disturbance 
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Is  determined  by  the  pilot  gain:  as  pilot  gain  is  increased,  the  cross¬ 
over  frequency  of  the  forward  transfer  function  is  raised.  Since  the 
parameters  of  Eq  ^8)  are  affected  by  the  plant  dynamics  and  the  distur¬ 
bance  input  in  the  region  of  crossover,  the  movement  of  the  crossover 
frequency  by  changes  in  pilot  gain  results  in  alternation  of  the  distur¬ 
bance  input  between  vg  and  pg  shown  in  Fig.  7.  In  this  case,  the  bandpass 
of  vg,  .38  radians/second,  and  the  bandpass  of  \>g,  .27  radians/second,  both 
give  an  effective  time  delay  of  .32  seconds  and  the  choice  of  a  predominant 
gust  is  probably  a  moot  point.  However,  in  determining  the  remaining  param 
eters  of  Eq  (8),  the  predominant  gust  input  will  be  identified. 

The  next  parameter  of  the  analytic  pilot  model  is  dependent  upon  the 
characteristics  of  the  plant.  Using  Cramer's  rule  again  on  Eq  (20),  the 
following  transfer  function  is  obtained: 


£  2.54(s  +  .78) (s  +  1.25) 

SPA  (s  *  ,19)(s2  +  1.77s  +1.14) 


(24) 


Figure  8,  a  log-magnitude  plot  of  this  transfer  function,  shows  that 
there  is  a  stretch  of  -20  decibel/decade  slope  from  a  frequency  of  1.25 
radians/second  to  a  frequency  of  7.25  radians/second.  This  information 
dictates  the  addition  of  a  pilot  lead  term  at  7.25  radians/second.  (Subse¬ 
quent  information  revealed  that  this  lead  term  was  not  necessary  for  stability 
as  was  originally  thought.  However,  the  term  was  carried  through  the  calcu¬ 
lations  and  will  be  included  here.)  At  this  point,  two.  of  the  three  required 
parameters  have  been  determined  and  Eq  (8)  has  the  form: 


Yp  =  Kpe"  • 32  ju,(l  +  s/7.25) 

Using  a  Fade  approximation  for  the  time  delay  term, 

KpCl  -  (.32/2)s](l  +  s/7.2) 
Cl + (.32/2)s] 


(25) 


(26) 
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Fig.  8.  Controlled  Element  Characteristics 
(Roll  Loop) 


6GC/EE/72-3 


/ 


Kp(l  -  .16s)(l  +  .14s) 
(1  +  .16s) 


.14Kp(s  +  7.25)(s-6.25) 


(s  +  6.25) 


(27) 


(28) 


Kp,  the  pilot  gain,  is  the  only  term  remaining  to  be  evaluated. 
McRuer  states  that  the  pilot  will  use  gain  to  obtain  a  crossover  fre¬ 
quency  for  the  forward  transfer  function  of  4.75  radians/second.  To 
determine  the  appropriate  gain,  the  following  requirement  holds: 

|YpYc|  =1  at  ui  =  4.75  radians/second  (29) 


From  Eqs  (24)  and  (25), 


(11.32) (1+1. 29s) (1+. 80s)  (l+.13s+. 01s2) (Kp) (l-.16s) (l+.14s) 
(1+5.223)(1+1.56s+.88s2)(1+.22s+.02s2)(1+.13s+.01s2)(1+.16s) 


=1 

s=juc 


(30) 


Evaluation  of  this  equation  reveals  that  at  a  frequency  of  4.75  radians/ 
second,  Kp  =  1.8  and  the  phase  angle  is  -205°.  Thus,  the  system  is  unstable 
and  the  pilot  must  operate  at  a  lower  frequency. 

Since  the  pilot  cannot  operate  at  the  desired  crossover  frequency,  he 
would  reduce  his  gain  to  achieve  a  crossover  frequency  with  an  acceptable 
stability  margin.  Thus,  any  pilot  gain  which  resulted  in  stable  operation 
would  be  acceptable. 

A  more  appealing  approach,  however,  is  to  choose  a  pilot  gain  which 
minimizes  roll  error,  the  pilot's  primary  concern  at  this  point.  The  pilot 
is  attempting  to  operate  the  closed  loop  system  in  Fig.  9,  where 
♦t  =  target  bank  angle  perturbation  (assumed  =  0), 

<t>  =  attacking  aircraft  bank  angle, 

Yp  =  pilot  model  given  by  Eq  (28),  and 
Yc  =  aircraft  dynamics  given  by  Eq  (24). 
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Fig.  9.  Roll  Angle  Control  System 


6PA  =  V  (31) 

=  -  Yp;/s  *  (32) 

Substituting  for  Yp  from  Eq  (28), 

.14K  (s-6.25)(s  +  7.19)p 

6PA=  —2 -  (33) 

(s)(s  +  6.25) 

Substituting  this  control  policy  Into  Eq  (20)  yields  the  following  matrix 
equation: 
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where  terms  A  through  E  are  given  In  Appendix  A,  and 

•  14K_(s -6.25XS +7.14) 

F  =  - 2 -  (35) 

(s)(s  +  6.25) 

Equation  (34)  represents  the  closed  loop  dynamics  of  the  control  system 
given  In  Fig.  9.  The  optimal  pilot  gain,  Kp  =  .7,  was  determined  by  vary¬ 
ing  Kp  In  Eq  (34)  from  .1  to  1.0,  and  then,  for  each  pilot  gain,  determining 
Eqs  (38)  and  (39)  and  evaluating  Eqs  (42)  and  ^43).  Figure  10  Is  a  graph 
of  the  data  for  pg  and  vg.  This  plot  presents  lata  for  Kp  vs  K*,  and  Kp  vs 
Kv  where 


"tv  =  Vv 

(36) 

“♦w  =  Vw 

(37) 

Using  Cramer's  rule,  the  following  transfer  functions  are  obtained  using 
a  pilot,  gain  of  .7  (Kp  =  .7  in  Eq  (35)): 


4  .00251(s)(s+.15)(s+.50)(s+3.84)(s-3.88) 
ni  (s+.27)2(s+.64)(s+1.74)(s2+2.47s+9.55) 


(38) 


J_  £  002527 ( s+. 80 )( s+l. 21) 

n2  (s+.64)(s+1.74)(s2+2.47s+9.55) 


(39) 


These  equations  relate  the  roll  angle  output  to  the  two  disturbance 
inputs;  since  the  predominant  disturbance  has  not  been  identified,  the 
effect  of  both  gust  inputs  must  be  evaluated. 

The  pilot  adjusts  his  gain  to  minimize  the  roll  angle  perturbations 
resulting  from  the  turbulence  inputs.  Since  the  gust  inputs  are  repre¬ 
sented  by  two  uncorrelated  white  noise  signals  of  unity  spectrum,  the 
spectral  densities  of  <p/m  and  <p/n2  are 


pg  EFFECT  (K*) 
vg  EFFECT  (Ky) 


!<„  (deg-sec/ft  *  104) 


o  p  p  o 

N  IO  U) 


(c0l  x  W39S-69P)  'bl 


PILOT  GAIN 
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9  (u)  =  4<M  .  ikM 

ni  rn 

N  V“)=^r“r 


The  root  mean  square  values  of  these  perturbation  power  spectra  are 


(«> 

(41) 


°+v  =^P'i,n,<“)  d“l  = 

X  4j[V2->  d"|  = 


=  .00428a,,  for  Eq  (36) 


=  .000570  ow  for  Eq  (37) 


(42) 


(43) 


Integral  (43)  was  evaluated  using  a  standard  fifth  order  integral 
equation  (Ref  21:4-73);  integral  (42)  was  not  of  the  proper  form  and  it 
was  evaluated  using  a  method  based  on  Simpson's  rule  (Ref  5:2-114). 

A  log-magnitude  plot  (Fig.  11)  of  the  forward  transfer  function 
reveals  that  for  Kp  =  .7,  the  crossover  frequency  is  2.0  radians/second 
with  a  phase  angle  of  -135°.  Thus  the  .system  is  stable;  the  pilot  has  a 
phase  margin  of  45°.  Further,  Fig.  8,  page  31,  shows  that  at  a  frequency 
of  2.0  radians/second,  the  plant  has  a  slope  of  -20  decibel/decade  and, 
consequently,  Eq  (9)  is  satisfied.  Also,  at  a  frequency  of  2.0  radians/ 
second,  Fig.  7,  page  29,  reveals  that  vg  is  the  primary  disturbance;  this 
is  confirmed  by  Fig.  10,  page  35,  which  shows  Ky  to  be  an  order  of  magni¬ 
tude  greater  than  Kw.  Finally,  the  complete  roll  pilot  model  is 

Y  .  -.0973(s  +  7.25)(s  -6.25)  ^ 

P  (s  +  6.25) 


The  development  of  this  pilot  model  has  been  explained  in  great 


detail  and  illustrates  the  procedure  which  will  be  applied  in  the  lateral 
and  longitudinal  compensatory  tracking  tasks.  Therefore,  the  parameters 
for  these  two  remaining  pilot  models  will  not  be  derived  in  detail. 


Fig.  11.  Forward  Transfer  Function  Characteristics 
(Roll  Loop) 
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Lateral  Tracking  v 

The  pilot  has  now  matched  the  roll  angle  of  the  target  and  thus  has 
Insured  a  stable  lead-pursuit  trajectory.  The  gust  disturbances,  however. 
Introduce  yaw  rate  perturbations  which  also  cause  the  plpper  to  be  dis¬ 
placed  laterally  from  the  target.  Therefore,  tiie  pilot  generates  another 
aileron  command  to  zero  the  lateral  tracking  error  u.  Again,  this  Is  a 
compensatory  tracking  task  and  an  analytic  pilo-;  model  must  be  developed 
for  this  case. 

The  pilot  observes  the  following  perturbation  errors: 

u  _  -  (.03516)$  (45) 

ni  "  nj  v  ; 


JL  -  ^  ~  (»03516H 
n2  H2 


(46) 


where  <|>/m ,  $/ni,  <J>/n2>  4>/n2  are  determined  from  Eq  (34)  using  Cramer's 
rule,  and 

y  =  $  -  (.03516)$  (47) 

Is  developed  by  substituting  z  and  a0  determined  in  Chapter  IV,  page  69, 
into  Eq  (7),  Chapter  II. 

Substituting  the  appropriate  values, 

JL  i  -07948(s  +  .15)(s  +  .19)(s  +  .5)(s2  +  2.98s  +  jJ2) 
ni  (s  +  .27)Z(s  +  .64)(s  +1.74)(s2  +  2.47s  +  9.55) 


y  _  _ .01325(s-  .  14 ) ( s  +  .52) _ 

nj  =  (s)(s  +  ,64)(s  +  1.74)(s2  +  2.47s  +  9.55) 

A  Bode  plot  of  these  two  equations.  Fig.  12,  reveals  that  vg  is  the 
dominant  gust  with  a  bandwidth  of  8.0  radians/second.  This  requires  that 
the  effective  time  delay  be  approximately  .12  seconds. 

The  plant  operating  characteristics  are 
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y  _  -  .03516$ 

6PA  6pa 


(50) 


Determining  the  appropriate  transfer  function*  from  Eq  (34)  and  substi¬ 
tuting  into  Eq  (50)  yields 


U  £  -  .2883(s  -  . 33 ) ( s  +  ,47)(;s  +  6.25) 

6PA  (s)(s  +  .64Xs  +  1.74)(s2  +  2.47s  +  9.55) 


(51) 


A  Bode  plot  of  this  transfer  function  shows  a  -20  decibel /decade 
nlope  from  u  =  0  to  to  *  .32.  Operation  Nove  this  range  of  frequency 
is  prohibited  by  the  stability  limits  of  the  plant. 

The  disturbance  input  bandwidth  and  the  plant  dynamics  given  by 
Eqs  (48)  and  (51),  respectively,  dictate  the  following  pilot  model  form: 

-Kp(s- 16.66) 


V 


(s  + 16.66) 


(52) 


As  In  the  preceding  case,  Kp,  the  pilot  gain,  is  selected  to  mini 
mize  closed  loop  tracking  error.  The#tracking  loop  model  is  shown  in 
Fig.  13. 

The  control  law  is 


where 


/+Kp(s-16.6)\  A  -  (Z-a0)f 


PA 


(s  +  16.6) 


=  Qip  -  Rj> 


lC<s-16.6) 

Q  =  - 

(s)  (s  +  16.6) 


K-Ts  -  16.6) 

R  =  — H -  (I-o0) 

(s)(s  +  16.6) 


(53) 

(54) 


(55) 


Substituting  Eq  (53)  into  Eq  (34),  the  following  matrix  equation  is 
generated: 
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Pilot 
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Disturbance  input 


Aircraft 
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Fig.  13.  Lateral  Tracking  Control  System 


s-Lp+ABL^-AFL^+ARL^  -Lp-CDLg^QALg^  -Lg-CEL^  p 

-Np+ABIi  ga~AFN  ga+ARN  ga  s-Nr-CDN^-QAN^  -N0-CEN6r  r 

■sYp-Y^  +ABYgI?-AFY6as+ARY5as  sYr-Y^CDY6rs-QAY6as  s2-sY6-CEY6rs  6 


-Lg-CEL6r 


CL6r 

AL6a 

Lp 

Lja 

CN6r 

AN6a 

Np 

Nr 

CY6rs 

AY*aS 

sYP 

sYr 

where  terms  A  through  E  and  the  appropriate  stability  derivatives  are 
given  in  Appendix  A,  term  F  is  given  in  Eq  (35),  and  terms  Q  and  R  are  given 
by  Eqs  (54)  and  (55),  respectively. 

As  in  the  roll  angle  control  task,  Cramer's  rule  is  applied  to  Eq  (56) 


\ 
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to  determine  the  closed  loop  determinant  and  the  transfer  functions 
relating  yaw  rate  and  roll  rate  outputs  to  the  disturbance  input.  From 
these  relationships,  an  equation  relating  y,  the  tracking  variable,  to  the 
disturbance  input  is  developed.  The  pilot  gain,  Kp,  is  then  iterated  from 
.1  to  1.0  in  this  equation  to  determine  the  gain  yielding  minimum  tracking 
error.  This  process  reveals  that  Kp  =  .1  minimizes  the  rms  intensity  of 
the -tracking  response  to  a  disturbance  input. 

Now,  at  a  pilot  gain  Kp =  .1,  the  pilot  model  is 

Y  _  - • l(s  -  16.66)  (57) 

P  (s  + 16.66) 

The  forward  transfer  function  of  the  system  depicted  in  Fig.  13  is 
obtained  by  combining  Eqs  (51)  and  (57): 


_y_  i  .0297(3  -  . 32) (s  +  ,47)(s  +  6.25) 
e  (s)(s  +  .64)(s  +  1.74)(s  +  2.47s  +  9.55) 


A  Bode  plot  of  Eq  (58)  reveals  that  the  crossover  frequency  is  .0004 

0 

radians/second.  Thus,  the  pilot  is  operating  on  a  -20  decibel/decade 
slope  of  the  plant  and,  consequently,  the  analytic  pilot  model  is  valid 
as  it  meets  all  constraints. 

The  closed  loop  system  depicted  in  Fig.  13,  and  mathematically  des¬ 
cribed  by  Eq  (56),  including  the  pilot  model  given  by  Eq  (57),  has  the 
following  response  characteristics: 

_y_  £  — . 00306( s  +  .15)(s  +  .19) (s  ♦  .5)(s2  ♦  2.99s  +9.09)  ^ 

ni  (s  +  .27)2(s  +  .64)(s  +  1.74)(s2  +  2. 48s  +  9.51) 

i>  .  -.00297(s +  .5)(s +  .15)(s  +  .18)(s2  +3.11s +9.84) 

_.  =  - - - ; -  (60) 

ni  (s+.27r(s+  .64)(s  +  1.74)(sz  +  2.48s  +9.51) 

$  .  +.00252(s-.001)(s  +  .i:)(s  +  .5)(s  +  3.81)(s  -  3.88) 

_  =  - - - 5 - '61/ 

ni  (s  +  .27)^(s  +  .64)(s  +  1.74)(s+ 2.48s +  9.51) 
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J_  i  +.00357(s  +  .15) (s2  +  .77s  +  .16)(s2  +  3.1s  +  9.86)  ^ 

ni  (s+  .27)2(s  +  .64)(s  + 1.74)(s2  +  2.48s  +9.51) 

These  four  equations  describe  the  propagation  of  disturbance  inputs 
through  the  system  dynamics  to  the  aircraft  controlled  variables.  Thus, 
they  are  powerful  analytic  tools  in  that  they  allow  us  to  quantitively 
evaluate  variations  in  ji,,  tbavtrao4i4ng •  variable,  and  <j>,  <p,  and  B,  the 
aircraft  controlled  variables  caused  by  a  specific  gust  input. 

Since  the  disturbance  is  a  random  variable  in  this  study,  it  must  be 
dealt  with  in  a  statistical  sense.  The  approach,  as  demonstrated  in  deter¬ 
mining  the  pilot  gains  in  the  two  analytical  pilot  models,  is  to  integrate 
the  power  spectrum  densities  of  the  filters  represented  by  Eqs  (59)  through 
(62)  and  then  to  relate  the  variance  of  the  output  to  the  variance  of  the 
input. 

Accomplishing  this  for  the  example  firing  situation  yields 


°w  = 

.1406 

°v 

(63) 

a*  = 

.2452 

av 

(64) 

- 

.1429 

av 

(66) 

aB  = 

.1759 

°v 

(66) 

Next,  the  magnitude  of  av  needs  to  be  determined.  From  Chalk  (Ref 
4:440-443)  the  expected  value  of  rms  vertical  gust  intensity  is 

ow  =  3.252  V(Pj)  (67) 

where 

Pj  =  probability  of  encountering  turbulence  at  15,000  ft 
=  .09  (from  Fig.  3,  page  443,  Ref  4) 

aw  =  .975  ft/sec  (68) 
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Also,  Eq  (23)  relates  av  and  ow;  thus, 

ov  =  1.39  ft/sec  (69) 

Substituting  Eq  (69)  into  Eqs  (63)  through  (66)  yields 

ctu  =  .1960°  (70) 

a*  =  .3418°  (71) 

a*  =  .1992°  (72) 

a0  =  .2452°  (73) 

One  final  task  must  be  accomplished  to  complete  the  model.  Fig.  4, 
page  22,  indicates  that  remnant  must  be  added  to  the  linear  element  of 
the  pilot  response  to  completely  duplicate  the  pilot's  performance. 
Remnant,  which  is  defined  as  the  portion  of  the  pilot's  control  output 
power  which  is  not  linearly  correlated  with  the  system  input,  contributes 
a  -relatively  large  portion  of  the  pilot's  response  in  a  task  such  as  the 
example  firing  situation  (Ref  16:131).  .  Levison  has  successfully  repre¬ 

sented  remnant  as  an  equivalent  observation  noise  source  (Ref  9).  For  a 
scalar  observation,  his  model  has  the  simple  structure: 


o*  =  P2(o2  +  a2) 
R  m 


where 


aR  =  variance  of  remnant  injected  entirely  as  an  equivalent 
visual  noise  source, 

P  =  constant  noise  ratio  independent  of  system  dynamics  and 
input  signal  bandwidth, 

a2  =  variance  of  displayed  signal,  and 
2 

a  =  variance  of  a  noise  source  accounting  for  the  minimum  or 
threshold  level  of  remnant. 

Thus,  the  lateral  tracking  error,  p,  should  be  revised  to  account  for 


the  effect  of  pilot  remnant. 


2  2.2 
°p  '  °4I+^R 
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where 

Oyj  -  Input  correlated  variance,  and 
o2R  "  remnant  Induced  variance. 

oj  =  o2n  +  P2(c2  +  a2 )  (76) 

=  (1  +  P2)oJ  +  P2o2  (77) 

For  the  F-4E  (R.?f  16:132), 

P  =  1.72 
om  =  .065° 

Substituting  Eqs  (70)  through  (73)  into  Eq  (77)  yields 


=  v'[(1  +  (1.72)2)(.19598)2  +  (1.72)2(  .065)2]  (78) 

=  .4057°  (79) 

=  .6894°  (80) 

=  .4119°  (81) 

aB  =  .5007°  (82) 


At  this  point,  a  model  relating  disturbance  input  to  the  output  of 
the  aircraft  controlled  variables  has  been  developed  and  applied  to  the 
example  firing  situation.  This  represents  accomplishment  of  the  first 
objective  of  this  chapter;  now,  a  similar  model  must  be  developed  for  the 
longitudinal  case. 

\ 

Longitudinal  Input  and  Controlled  Element 

The  vertical  turbulence  spectrum,  which  induced  a  roll  rate  in  the 
lateral  case,  also  disturbs  the  aircraft  aerodynamical ly  in  the  longitu- 
dinal  axis.  This  vertical  gust,  w  ,  produces  angle  of  attack,  pitch  rate, 

o 

and  angle  of  attack  rate  perturbations  which  are  observed  by  the  pilot  as 
disturbances  in  the  pitch  attitude  of  the  aircraft.  Thus,  when  the  attacking 
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aircraft  is  perturbed  from  its  nominal  lead-pursuit  trajectory,  the  pilot 
observes  vertical  displacement  of  the  plpper  from  the  target  aircraft  as 
well  as  the  lateral  ploper  displacement  discussed  in  the  previous  section. 

As  in  the  lateral  case,  the  vertical  gust,  wg,  can  be  modeled  by 
passing  a  Gaussian  random  "white"  noise  signal  through  a  low-pass  filter 
(Ref  4:460-461)  The  filter  characteristics  can  be  modeled  as  shown  in 
Fig.  14,  where 

v'OwMJnlCl  +  V(3)L*s/U0]ow 

TWg  =  - 2 -  (83) 

(1  +  LmS/Uq) 

n3  s  white  noise  signal , 
s/Uq 

T0P  = - 

8  [l  +  (4b/rrU0)s]  (84) 

U0,  L^,  and  b  are  given  after  Eq  (3). 

Substituting  the  appropriate  values  for  the  example  firing  situation 
into  Eqs  (83)  and  (84)  yields 

a_  -.0007572(s  +  .31)ou 

A  =  - - - -  (85) 

n3  (S  +  .54)2 


ap  -.01473(s2  +  .31s)aw 

“  =  - i -  (86) 

n3  (S  + .54)Z(s +19.46) 

.01473(s2  +  .31s)o 

^  =  - — - -  (»7) 

n3  (S+  .54)  (s  +  19.46) 

The  approximate  equations  of  aircraft  longitudinal  motion,  the  "short 
period"  equations,  describe  a  portion  of  the  controlled  element  as  well  as 
additional  filtering  for  the  gust  input.  Hence,  these  equations  relate 
the  aircraft  longitudinal  response  to  elevator  and  vertical  gust  inputs 
(Ref  1:32): 


Fig.  14.  Model  of  the  Effect  of  Vertical  Gusts  on  the  Longitudinal 
Dynamics  of  an  Aircraft  (Adapted  from  Ref  17:76) 


s-Mg  -t%s-Ma 

0 

Ms  Me  Ma  Mi 

-l-Zg  S-Za 

La 

^6  2g  z0  0 

“g 

“g 


(81 


The  aircraft  dynamics  are  altered  by  the  stability  augmentation  sys 
tern  (SAS)  incorporated  in  the  pitch  channel  and  hence  the  effects  of  the 
pitch  SAS  must  be  included  in  Eq  (88). 

The  pitch  axis  control  system  can  be  represented  as  (Ref  17:155) 

4(s)  =  Gp(s)[K<s)5pE(s)  +  Hp(s)0(s)J  (8! 
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6(s)  =  elevator  deflection, 

Gp(s)  =  power  cylinder  dynamics, 

K(s)  =  linkage  dynamics, 

Ap^Cs)  =  longitudinal  control  stick  deflection,  and 
Hp(s)  =  pitch  rate  feedback  dynamics. 


Substituting  Eq  (39)  into  Eq  (88)  yields 


s-M$-M6GH  -MaS-Ma 

0 

M6GK 

Me  Ma  Md 

— 1— Z^— Z^GH  s— Zo 

Act 

z6gk 

ze  za 

“g 


(90) 


The  above  matrix  equation  is  a  mathematical  model  of  the  longitudinal 
axis  of  the  attacking  aircraft  which  the  pilot  must  control  to  success¬ 
fully  track  the  target.  Cramer's  rule  can  be  applied  to  Eq  (90)  to  obtain 
transfer  functions  relating  the  aircraft  controlled  variables,  e  and  Aa, 
to  both  gust  and  pilot  inputs.  These  equations  can  be  used  to  quantitively 
evaluate  the  perturbati  >ns  in  0  and  Aa  caused  by  the  turbulence  input  and 
also  to  determine  the  pilot's  actions  necessary  to  restore  the  aircraft 
to  the  nominal  lead-pursuit  trajectory.  The  numerical  entries  for  each 
element  of  this  matrix  equation  are  listed  in  Appendix  A  and  will  be  sub¬ 
stituted  in  the  appropriate  positions  in  the  succeeding  equations. 


Longitudinal  Tracking 

The  attacking  aircraft  has  been  perturbed  from  its  nominal  lead-pursuit 
trajectory  by  gust  inputs;  the  pilot  observes  this  perturbation  as  both 
lateral  and  longitudinal  displacement  of  the  pipper  from  the  target.  In 
the  preceding  sections,  a  model  has  been  developed  to  describe  pilot  control 
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of  lateral  pipper  displacement;  a  similar  model  will  now  be  developed  to 
account  for  control  of  longitudinal  pipper  replacement  by  the  pilot. 

As  discussed  in  Chapter  II,  the  pilot  controls  &,  the  angle  in  the 
pitch  axis  between  the  pipper  and  the  target,  which  is  generated  by  the 
gust  Input.  This,  the  pilot  is  confronted  with  a  third  compensatory 
tracking  task.  This  requires  the  development  of  another  analytical  pilot 
model,  substitution  of  this  model  into  the  system,  and  a  loop  closure  to 
control  the  pitch  tracking  variable  0. 

The  closed  loop  system  describing  the  longitudinal  tracking  task  is 
shown  in  Fig.  15. 

Before  proceeding,  it  is  necessary  to  determine  the  transfer  func¬ 
tion  relating  e,  the  aircraft  pitch  rate,  to  0p,  the  Dipper  position 
on  the  combining  glass.  The  equations  describing  the  dynamics  of  the 
fire  control  system,  including  the  sight,  are  (Ref  15:17-19) 

9P  _  as2  +  bs  +  c  (gij 

0  s(s  +  AV/DHTfS  +  1) 

where 


e  =  pitch  rate  of  the  aircraft, 

9p  =  depression  of  pipper  below  zero  depression  line, 

AV  =  VT  -  VA  (92) 

with  VA  =  attacking  aircraft  velocity 
VT  =  target  aircraft  velocity 
D  =  target  present  range, 

Tf  =  time  of  flight  of  projectile, 

a  =  (1  +  WTf)  (93) 

b  =  (AV/DXl  -  WTf  +  ZttTf)  -  Z0(l  +  VTTf/D)  (94) 
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Fig.  15.  Longitudinal  Tracking  System 


c  =  -  ZaVT/D  *  (95) 

W  =  -  (VA/DVf)(Za  +  2JV/Tf)  (96) 

The  aerodynamic  constants  required  in  these  equations  are  found  in 
Appendix  A;  the  firing  situation  parameters  and  ballistic  solution  data 
are  found  in  Chapter  IV.  Substituting  the  appropriate  values  results 
in  the  following  equation: 

6_P  _  .825(s  +  ,47)(s  +  1.65)  ^ 

®  (s)(s  +  1.32) (s  -  .053) 

Applying  Cramer's  rule  to  Eq  (90)  to  determine  the  transfer  function 
e/n3  and  combining  this  equation  and  Eq  (97)  yields  Eq  (98),  which  describes 
the  perturbation  observed  by  the  pilot: 


*In  the  referenced  technical  report,  Eq  (95)  is  preceded  by  a  positive, 
not  by  a  negative  sign  as  in  this  report.  The  sign  change  was  made  for 
two  reasons.  A  negative  sign  coupled  with  a  negative  Za  results  in  a 
numerically  positive  term  which  is  required  for  stability.  Second,  when 
the  contractor's  data  was  substituted  into  his  equations,  a  numerically 
positive  term  occurs  in  his  following  equations.  This  can  only  occur  if 
Za  is  positive  or  if  both  the  term  and  Za  are  negative. 
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i  .0099S5(s  +  .31)(s +  .47)(s  +  1.65) _ 

ns  (s-  .053)(s  +  .54)2(s +1.32)(s2  + 7.28s + 18.10) 


(98) 


A  Bode  plot  of  this  transfer  function  shows  the  input  bandwidth  to  be 
.05  radian/second. 

The  plant  transfer  function  is 


irr- 

PE 


19. 8422(s  +  .47) (s +  ,98)(s  +  1.65) 

(s)  (s  —  .053)(s  +1.32)(s  +7.:  8s +  18.10) 


(99) 


A  Bode  plot  of  Eq  (99)  indicates  that  the  pilot  would  add  a  lag  at  1 
radian/second  in  order  to  obtain  a  -20  decibel/decade  slope  from  .47 
rad/sec  to  1.32  rad/sec  and  from  1.65  rad/sec  to  4.25  rad/sec. 

From  the  information  contained  in  Eqs  (98)  and  (99),  Eq  (8)  becomes 

-  K_.(s  -  6.25) 

Y  = - 2 -  (100) 

(s +l)(s  +  6.25) 


The  pilot  will  select  his  gain  to  minimize  pitch  tracking  error. 

As  in  the  lateral  case,  the  optimum  pilot  gain,  Kp,  is  selected  by  closing 
the  tracking  loop  and  iterating  Kp  to  determine  the  minimum  root  mean 
square  tracking  error. 

Thus  the  control  law  for  the  system  depicted  in  Fig.  15  is 

.825K  (s  -  6.25)<s  +  .47)(s  +  1.65)9 

6  =  - £ -  (101) 

(s)(s  +  1.32)(s  -  .053)(s  +  l)(s  +  6.25) 

Iterating  pilot  gain  in  Eq  (101)  from  .1  to  1.0,  the  optimal  gain  is 
determined  to  be  Kp  =  .4.  Substituting  Eq  (101)  into  Eq  (90)  and  applying 
Cramer's  rule,  we  have,  for  Kp  =  .4, 

fP  f  .03725(s) (s  -  ,05)(s  +  .31) (s  +  l)(s  +  1.32)  ^1Q2j 

n3  (s+  .54)2(s2  +  .23s  +  .40)(s2  +2. 02s  +  1.03)(s  +  2.91) 

For  this  gain,  the  rms  tracking  error  is 
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00  =  .08855ow  =  .0863° 


(103) 


Adding  pilot  remnant, 


o0  =  .1850° 


(104) 


To  determine  the  operating  frequency  of  the  closed  loop  system,  the 
crossover  frequency  of  the  forward  transfer  function  must  be  determined. 
The, .forward  transfer  function,  with  pilot  gain,  Kp,  unspecified,  is  a 
combination  of  Eqs  (99)  and  (100), 


_ 19.8422(s  +  .47)(s  +1.65)(s-6.25)(Kp) 

(s)(s-  .053)(s  +1.32)<s  +  7.28s  +  18. 10)(s  +  6.25) 


(105) 


A  Bode  plot  of  Eq  (105),  Fig.  16,  with  Kp  =  .h,  reveals  that  the  pilot 
Is  operating  where  the  system  slope  is  -40  decibels/decade.  This  is  an 
unacceptable  operating  point,  even  though  tracking  error  is  at  a  minimum. 
Equation  (99)  indicated  that  a  slope  of  -20  decibels/decade  exists  in 
the  frequency  range  of  .47  to  1.32  radaans/second.  The  Bode  plot  of  Eq 
(105)  reveals  that  a  pilot  gain  Kp=  .8  will  give  a  crossover  frequency 
of  .90  radians/: econd.  This  is  exactly  in  the  middle  of  the  -20  decibels/ 
decade  slope  and  hence  acceptable  to  the  pilot  as  an  operating  point. 

For  Kp  =  .8,  Eq  (102)  becomes 

!£  =  -001668(s)(s  -  .05)  (s  +  .31)(s  +  1)  (  s  +  1.32)  /iQ6j 

n3  (s  +  .54)  (s  +  .72)(s  +  .95)(s  +  2.02)(s  +.23s+1.22) 


Integration  of  this  equation  and  addition  of  pilot  remnant  results 


in  the  following  root  mean  square  tracking  error: 


Similarly, 


oe  =  .2482° 
aLa  =  .2112° 


(107) 

(108) 


Thus,  in  securing  an  acceptable  operating  regime,  mean  square  tracking 
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error  has  risen  considerably  from  Eq  (104)  to  Eq  (107). 

Equation  (106)  represents  completion  of  the  second  model;  It  is  a 
mathematical  statement  which  relates  pitch  tracking  error,  e,  to  the  gust 
Input  when  the  pilot  is  attempting  to  return  the  pipper  to  the  target  and 
hence  reduce  tra;king  error  to  zero. 

To  this  point,  the  lateral  case  and  the  longitudinal  case  have  been 
dealt  with  as  separate  entities.  However,  only  one  pilot  controls  both 
of  these  channels  and,  hence,  the  interdependence  of  these  two  tasks 
must  be  considered. 

Homogeneous  Pilot  Model 

When  the  pilot  is  confronted  with  a  two-axis  tracking  task,  he  will 
use  the  same  control  policy  in  both  channels  (Ref  10).  Thus,  a  single 
pilot  model  must  be  developed  to  replace  the  lateral  pilot  model,  Eq  (57), 
and  the  longitudinal  pilot  model,  Eq  (100),  with  Kp  =  .8.  The  same  adjust¬ 
ment  rules  are  used  to  determine  the  parameters  of  this  homogeneous  pilot 
model;  additionally,  the  model  is  constrained  to  cause  the  system  cross¬ 
over  frequency  to  occur  along  a  -20  decibel/decade  slope  in  both  axes. 

The  bandwidth  of  vg  predominates,  hence  the  lateral  effective  time 
delay  will  be  used  in  the  homogeneous  model.  The  lag  at  one  second  in 
the  longitudinal  pilot  model  will  not  be  included  in  the  homogeneous  model 
because  it  greatly  increases  lateral  tracking  error.  Further,  close  inves¬ 
tigation  of  the  longitudinal  case  reveals  that  the  system  goes  unstable 
at  Kp  =  1 . 1  with  the  single  axis  pilot  model.  Thus,  the  pilot  is  unable 
to  take  advantage  of  the  extended  -20  decibel/decade  slope.  Removing  the 
lag  would  cause  the  pilot  to  operate  in  a  somewhat  more  confined  -20 
decibel/decade  region;  a  benefit,  however,  is  that  longitudinal  tracking 
error  is  reduced  when  the  lag  is  omitted.  Finally,  in  order  to  operate 
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In  the  middle  of  a  -20  decibel/decade  region  in  the  longitudinal  case, 
a  pilot  gain  of  K.0=  .6  must  be  used;  at  this  gain,  the  lateral  system 
is  also  operating  on  the  desired  slope. 


The  above  constraints  yield  the  following  homogeneous  pilot  model 


-.6(s  -  16.667) 
(s  +  16.667) 


(109) 


'•Returning  to  the  preceding  section  on  longitudinal  tracking,  Eq 
(100)  is  replaced  by  Eq  (109)  and  the  computations  of  that  section  are 
repeated  with  the  fixed  gain  homogeneous  pilot  model.  The  root  mean 
square  values  of  the  longitudinal  controlled  variables,  with  pilot  rem¬ 
nant  included,  are  determined  to  be 


oq  =  .1308° 
=  *1443° 


(no) 

(111) 


In  addition  to  the  homogeneous  pilot  model,  another  important  effect 
must  be  considered  in  the  two-axis  tracking  task.  The  larger  lateral 
tracking  errors  interfere  with  the  pilot's  ability  to  sense  error  and 
error  rate  in  the  longitudinal  axis;  consequently,  longitudinal  tracking 
performance  is  degraded  from  that  predicted  by  single  axis  analysis  even 
with  homogeneous  pilot  models  and  remnant  included.  This  degradation  is 
called  the  visual-motor  interference  and  can  be  expressed  analytically  as 


(Ref  16:131-132) 


oeT  =  y/toQ  +  K^oJ] 


(112) 


wnere 


aeT  =  two  axis  longitudinal  tracking  error  intensity, 

cry  =  single  axis  lateral  tracking  error  intensity, 

Kfl  =  interference  coefficient  =  .558. 
u 
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Applying  this  relationship  to  Eq  (110)  and  Eq  (111)  yields 

ce  =  .2615°  (113) 

°Ac.t  =  '2686°  <114) 

Similarly,  substituting  Eq  (109)  instead  of  Eq  (57)  into  Eq  (56)  and 
recomputing  the  root  mean  square  values  of  the  lateral  controlled  variables 
and  the  lateral  tracking  error  yields 


- 

.4148° 

(115) 

a*  = 

.6824° 

(116) 

°6  = 

.5007° 

(117) 

= 

.4059° 

(118) 

Investigation  of  the  forward  transfer  functions  reveals  that  the 
operating  frequencies  are  .72  radians/second  in  the  longitudinal  case  and 
.0012  radians/second  in  the  lateral  case. 

The  six  equations  listed  above  are  the  completed  longitudinal  and 
lateral  tracking  error  models,  respectively,  for  the  example  firing  situa¬ 
tion.  Thus,  for  a  given  level  of  gust  inputs,  the  tracking  error  in  both 
axes  has  been  predicted.  The  final  task,  to  be  accomplished  in  the  next 
chapter,  is  to  quantitively  evaluate  the  terminal  miss  distance  caused  by 
these  tracking  errors. 

Special  Longitudinal  Case 

Occasionally,  the  pilot  is  forced  to  track  and  fire  with  a  fixed 
sight  instead  of  the  lead-computing  sight.  In  this  situation,  the  sight 
display  indicates  a  correct  firing  solution  for  only  one  case  as  the  pipper 
remains  depressed  at  a  constant  angle  regardless  of  the  state  of  either  the 
target  aircraft  or  the  attacking  aircraft.  For  the  purposes  of  this  study, 
it  will  be  assumed  that  the  fixed  depressed  reticle  is  indicating  the  correct 
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aim  point  in  the  example  firing  situation;  the  object  of  this  investiga¬ 
tion  is  to  remove  the  dynamics  of  the  fire  control  system  and  to  determine 
the  effect  that  this  action  will  have  on  the  longitudinal  and  lateral 
tracking  performance. 

The  analysis  proceeds  by  removing  the  fire  control  components  from 
the  system  illustrated  in  Fig.  15  and,  consequently,  omitting  the  terms 
of  Eq  (-97)  from  the  appropriate  entries  of  Eq  (90).  Paralleling  the  pro¬ 
cedure  contained  in  the  Longitudinal  Tracking  section,  calculations  reveal 
that  the  pilot  observes  a  disturbance  input  with  a  bandwidth  of  1.0  radians/ 
second  and  the  plant  has  a  -20  decibel/decade  slope  from  a  frequency  of 
0.0  radians/second  to  a  frequency  of  .98  radians/second. 

This  information,  coupled  with  the  information  derived  in  the  Lateral 


Tracking  section,  dictates  the  following  homogeneous  pilot  model: 

Y  _  -  .Ks-16.66) 
p  '  (s  + 16.66) 


(119) 


Repeating  the  computations  outlined  in  the  preceding  section  with 
the  pilot  model  given  by  Eq  (119),  the  following  root  mean  square  values 
are  predicted  for  longitudinal  and  lateral  tracking  error: 


o0  =  .2656° 


oy  =  .4057° 


(120) 

(121) 


Bode  plots  of  the  forward  transfer  functions  reveal  that  both  closed 
loop  systems  are  operating  on  -20  decibel/decade  slopes  at  frequencies  of 
.08  radians/second  in  the  longitudinal  case  and  .0004  radians/second  in 
the  lateral  case. 
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IV.  Mathematical  Model  of  the  Aerial  Gunnery  Task 


Introduction 

In  this  chapter  the  prediction  angle  will  be  defined  and  the  parameters 
of  the  nominal  lead-pursuit  trajectory  will  be  calculated  for  the  example 
firing  situation.  Then,  the  attacking  aircraft  will  be  disturbed  from 
its  nominal  state  by  gust  inputs.  The  cannon  will  be  fired  from  this 
perturbed  state  and  the  projectile  terminal  miss  distance  resulting  from 
the  errors  in  the  tracking  and  controlled  variables  of  the  perturbed  state 
will  be  determined. 


Prediction  Angle 

Description.  In  order  to  hit  a  moving  target  with  an  aircraft-mounted 
cannon,  it  is  necessary  to  create  the  proper  angle  between  the  tracking 
line  and  the  weapon  line.  This  angle  is  called  the  prediction  angle  and 
consists  of  the  following  components  (Ref  22:2) 

•  lead 

•  curvature  correction 

•  jump  correction 

Lead .  For  the  projectile  to  strike  the  target,  the  component  of  the 
projectile  velocity  perpendicular  to  the  tracking  line  must  be  equal  to 
the  component  of  the  target  velocity  perpendicular  to  this  same  line. 

Figure  17  presents  the  two  dimensional  case. 

From  Fig.  17, 


Vp  sin  L  =  V,p  sin  p 


Hence, 


L  =  arc  sin 


VT  sin  p 


(122) 


(123) 
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where 

L  =  lead  angle, 

Vp  =  projectile  total  velocity 
p  =  target  angle  off,  and 
VT  =  target  velocity  (Ref  2:3) 

Curvature  Correction.  There  are  three  components  of  curvature 
correction: 

The  wind  curvature  angle  is  zero  in  the  aerial  gunnery  situation  as 
the  attacking  aircraft,  the  target  aircraft,  and  the  projectile  are  all 
moving  in  the  same  air  mass.  In  firing  at  a  ground  target,  however,  a 
left  crosswind  would  cause  the  ground  track  of  the  projectile  to  move  to 
the  right  as  it  travels  to  the  target.  In  the  example  firing  situation, 
wind  curvature  is  ignored. 

The  drift  correction  angle  is  caused  by  the  interaction  of  the  torque 
produced  by  the  aerodynamic  lift  acting  on  the  projectile  ?  .d  the  angular 
momentum  of  the  projectile  due  to  spin.  The  result  of  this  action  is  a 
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"drift"  of  the  projectile  to  the  right  If  the  spin  is  to  the  right.  The 
effect  of  this  interaction  is  small  and  can  be  neglected  (Ref  2:6). 

The  gravity  drop  angle  is  the  result  of  the  action  of  gravitational 
acceleration  upon  the  projectile  during  its  time  of  flight.  The  situation 
is  depicted  in  Fig.  18  (Ref  2:7). 


vertical 


Actual 

trajectory 


"tracking  li°e 


Fig.  18. 

Gravity  Drop  Geometry 

From  Fig.  18, 

A 

■  1 
"  2 

B 

=  7  STf  cos  E 

Hence, 

sin  C 

1  2 

2  gt f  cos  E 

R 

and 

R 

tf 

1 

<1 
$  1 

Thus, 

C 

=  arc  sin  rR  cos  E 

2V2 

2  PA 
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where 


tf  =  projectile  time  of  flight, 
g  =  gravitational  acceleration  constant, 

E  =  elevation  angle  of  weapon  line  above  the  local  horizontal, 

R  =  target  range, 

C  =  gravity  drop  angle,  and 
VpA  =  average  projectile  velocity. 

Jump  Correction.  There  are  two  component:,  of  the  jump  correction 
angle: 

The  windage  jump  angle  is  created  by  the  misalignment  of  the  weapon 
line  and  the  aircraft  velocity  vector.  As  the  projectile  leaves  the 
muzzle  and  enters  the  airstream,  it  has  an  angle  of  attack  and  experiences 
a  net  lift  force*  This  lift  force  inclines  the  true  projectile  path  away 
from  the  weapon  line  towards  the  aircraft  velocity  vector.  This  term  can 
be  neglected  in  the  aerial  gunnery  situation  (Ref  2:3). 

The  velocity  jump  angle  is  a  consequence  of  the  fact  that  the  total 
projectile  velocity  is  the  vector  sum  of  the  muzzle  velocity  vector  and 
the  aircraft  velocity  vector.  Figure  19  is  a  geometric  description  of 
velocity  jump  (Ref  2:4). 

From  Fig.  19, 

smJ-.f--  S2SSJ.  (129) 

P  vp 


arc  sin 


sin  £  V 


(130) 


where 


J  =  jump  angle, 

£  =  angle  between  the  muzzle  velocity  vector  and  the  aircraft 
velocity  vector. 
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V  =  aircraft  velocity  vector,  and 
vp  =  projectile  total  velocity. 

Determination  of  the  Nominal  Trajectory 

Introduction.  The  parameters  of  the  nominal  lead-pursuit  trajectory 
are  determined  by  an  iterative  algorithm  on  the  digital  computer.  This 
method  of  solution  is  dictated  by  the  fact  that  no  closed-form  solution 
of  the  prediction  angle  equation  exists  when  projectile  drag  and  aircraft 
angle  of  attack  perturbations  are  included  in  the  mathematical  model,  and 
these  latter  effects  must  be  considered  if  the  example  firing  state  is  to 
be  a  realistic  situation. 

The  calculation  of  the  firing,  state  parameters  will  be  described 
briefly  in  this  section.  The  intent  is  to  follow  the  mathematical  flow 
of  the  digital  computer  program  explaining  in  detail  unusual  or  interesting 
aspects  of  the  computations.  Finally,  the  parameters  of  the  nominal  lead- 
pursuit  trajectory  will  be  stated. 
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Calculation  of  Prediction  Angle.  The  calculations  are  initiated  by 
establishing  a  tracking  line  from  the  attacking  aircraft  to  the  target 
aircraft.  Then  the  target  aircraft's  velocity  vector  is  resolved  into 
the  tracking  line  reference  system  described  in  Chapter  II.  Once  the 
velocity  components  of  the  target  aircraft  are  known  in  the  tracking  line 
reference  system,  Eq  (123)  is  applied  in  the  TLYand  TLZaxes  to  determine 
,.az-imuth  and. .eleve frion  ‘  lead  angles,  respectivel /. 

Then,  the  angle  of  attack  necessary  to  mantain  level,  unaccelerated 
flight  in  the  example  firing  situation  is  substituted  into  Eq  (130)  to 
determine  the  jump  correction  in  the  elevation  plane. 

There  are  two  noteworthy  aspects  of  the  jump  correction  angle. 

First,  the  angle  5  must  be  defined  and  measured  with  care.  As  illus¬ 
trated  in  Fig.  20,  two  distinct  situations  can  occur. 

In  the  first  case,  the  weapon  line  is  between  the  tracking  line  and 


_ CASE  1  _ | _ CASE  2 _ | 

Fig.  20.  Jump  Angle  Condition 

the  aircraft  velocity  vector  and  the  jump  angle  is  subtracted  from  the 
total  lead  angle;  in  the  second  case,  the  weapon  line  does  not  lie  between 
the  tracking  line  and  the  aircraft  velocity  vector  and  the  jump  angle  is 
added  to  the  total  lead  angle.  Consequently,  when  the  cannon  muzzle  is 
depressed  below  the  aircraft  velocity  vector,  t,  will  be  considered  negative 
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and  when  the  weapon  line  Is  elevated  above  the  aircraft  velocity  vector, 

5  will  be  considered  to  be  positive. 

Second,  it  is  assumed  that,  in  the  nominal  firing  situation,  the 
tracking  line,  the  weapon  line,  and  the  attacking  aircraft's  velocity 
vector  all  lie  ir  the  elevation  plane  of  the  attacking  aircraft;  thus, 
the  jump  angle  correction  occurs  in  the  aircrart  longitudinal  plane, 
how,  since  the  total  lead  angle  is  created  and  maintained  by  pitch  control, 
the  jump  correction  angle  is  either  added  to  or  subtracted  from  the  total 
lead  angle.  It  should  be  noted,  however,  that  the  attacking  aircraft's 
velocity  vector  can  be  displaced  out  of  the  longitudinal  plane  by  intro¬ 
ducing  sideslip.  The  F-4  fire  control  system  cannot  compensate  for  this 
condition.  Thus,  this  situation  is  an  error  producing  state  and  will  be 

considered  to  be  a  perturbation  from  the  nominal  firing  position. 

0 

Finally,  the  gravity  drop  is  calculated  using  Eq  (128).  Note  that 
in  the  first  iteration,  the  projectile  initial  velocity  rather  than  the 
projectile  average  velocity  is  used  in  Eq  (128). 

The  vector  addition  of  these  three  angles  forms  the  prediction  angle; 
however,  this  first  prediction  angle  must  be  corrected  for  projectile  drag 
and  aircraft  angle  of  attack. 

Correcting  the  Prediction  Angle  for  Projectile  Drag.  Aerodynamic 
drag  acts  to  continuously  reduce  the  velocity  of  the  projectile  during 
its  time-of-flight.  This  velocity  reduction  in  turn,  alters  the  lead 
angle  and  the  gravity  drop  terms.  Therefore,  it  is  necessary  to  deter¬ 
mine  the  projectile  average  velocity  and  time  of  flight  in  order  to 
calculate  the  correct  prediction  angle. 

The  drag  force  on  the  projectile  is  (Ref  17:149-150) 

CdpV2S 
D  =  — - 


(131) 
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where 

Cp  =  coefficient  of  drag, 

S  =  suitable  scaling  area, 
p  =  atmospheric  density,  and 
V  =  projectile  speed. 


From  Mechanics, 


F  =  ma 

(132) 

C^S 

ma  = - 

2 

(133) 

CjOV^S 
a  "  2m 

(134) 

where 

m  =  projectile  mass,  and 

u  -'protfecti te  acceleration. 

KpV2 

2 

(135) 

where 

<TjS 

k  =  -t- 

(136) 

a  =  deceleration  of  projectile. 

Now, 

.  dV  _  KpV2 

3  "  dt  2 

(137) 

Separating  variables. 

dV  Kp 

v*  T d 

(138) 

Integrating, 

1  _  Kpt  1 

"  V  "  "  2  "  Vo 

(139) 
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Hence, 


Also, 


Integrating, 


Rearranging, 


(2  +  KpV0t) 


V  -  <«  .  ;v° 

dt  2  +  KpVot 


‘(rr 


KpVot  ( 


2  ( 

R  =  —  In  1  +  - - 

Ko  \  2 


(140) 


(141) 


(142) 


(143) 


(144) 


t  _  2e(KpR/2)  -  2 


(145) 


V  =  Voe 


-fKpR/2) 


(146) 


At  this  point,  K  must  be  determined.  Rearranging  Eq  (144)  and  substituting 


for  the  following  variables  (Ref  14:49): 


where 


K  =  In \—j  =  .2011083  ft2/slug 


V0  =  4000  ft/sec, 

V  =  2000  ft/sec  at  R  =  2900  ft,  and 
p  =  .002377  slug/ft  at  sea  level. 


(147) 


To  correct  the  prediction  angle  for  projectile  drag,  the  bullet  time- 
of-flight  is  determined  from  Eq  (145)  and  the  final  velocity  from  Eq  (146). 
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An  average  velocity  is  then  determined  and  an  augmented  range  is  formed 
by  adding  the  target  distance  covered  during  the  projectile's  time-of- 
flight  to  the  original  distance.  This  iterative  procedure  is  continued 
until  two  successive  augmented  ranges  meet  some  stopping  criterion.  The 
average  velocity  is  then  used  to  create  another  lead  angle  and  the  new 
time-of-fl ight  and  average  velocity  are  used  to  compute  gravity  drop, 
ffu w,  since  the ‘increased  >eud  angle -has  reduced  the  bullet's  forward 
velocity,  another  iterative  loop  is  activated  which  corrects  the  predic¬ 
tion  angle  for  drag  again.  This  outer  loop  runs  until  two  successive 
times-of-f light  are  within  some  stopping  criterion.  At  this  point,  the 
prediction  angle  is  corrected  for  projectile  drag. 

Correcting  the  Prediction  Angle  for  Angle  of  Attack.  The  aircraft 
angle  of  attack  affects  £,  the  angle  between  the  weapon  line  and  the 
attacking  aircraft's  velocity  vector.  Consequently,  the  jump  angle  and 
the  total  lead  angle  must  be  corrected  for  angle  of  attack. 

Angle  of  attack  can  be  calculated  from  the  attacking  aircraft's 
state.  First,  an  aircraft  weight  is  selected;  then  the  lift  required  to 
maintain  equilibrium  is  determined: 

L  cos  9 


W  = 


cos  0 


where 


L  = 


CLjPV2S  . 


(148) 


(149) 


The  previously  generated  prediction  angle  yields  the  aircraft  bank 
angle,  $,  and  pitch  attitude,  o.  The  aircraft  speed  is  specified,  density 
is  determined  from  the  specified  altitude,  and  S  is  known  for  the  F-4E. 
This  allows  us  to  solve  Eq  (148)  for  C^,, 

2W  cos  0 


'bp 


cos  $  CV2p 


(150) 
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For  the  particular  attacking  aircraft  state,  curves  of  C^,  vs  a 
will  yield  the  angle  of  attack.  Since  the  cannon  is  depressed  3  degrees 
below  the  wing  zero  lift  line  (Ref  11), 

£  =  cx  -  3°  (151) 

The  value  of  s  is  used  to  generate  a  corrected  prediction  angle;  this 
new  prediction  angle  yields -a  different  angle  of  attack.  Consequently,  an 
iterative  procedure,  using  the  closeness  of  two  succeeding  values  of  £  as 
a  stopping  criterion,  is  used  to  correct  the  prediction  angle  for  attacking 
aircraft  angle  of  attack. 

The  method  used  to  determine  the  coefficient  of  lift  (Cl)  implies 
that  available  aircraft  thrust  is  sufficient  to  offset  all  lift  loss  due 
to  pitch  attitude.  This  is  true  only  for  a  limited  range  of  pitch  attitude, 
G- load,  airspeed  states.  Exceeding  this  envelope  will  result  in  aircraft 
deceleration. 

Further,  the  attacking  aircraft  eft id  the  target  aircraft  are  considered 
to  be  in  equilibrium.  This  state  is  valid  only  in  a  small  number  of  aerial 
gunnery  situations.  Nevertheless,  situations  do  occur  where  the  attacking 
aircraft  has  excess  thrust  and  the  target  is  in  unaccelerated  flight  and, 
hence,  this  situation  is  a  valid  example  firing  situation. 

Computed  Example  Firing  Situation.  Once  the  prediction  angle  is 
corrected  for  projectile  drag  and  aircraft  angle  of  attack,  the  parameters 
of  the  static  firing  situation  can  be  determined.  The  aircraft  geometry, 
angle  of  attack,  and  the  prediction  angle  can  be  used  to  determine  the 
aircraft  pitch  attitude  in  stability  axes.  The  roll  angle  control  require¬ 
ment  dictates  the  bank  angle.  The  prediction  angle  iterative  solution 
generates  an  entire  set  of  projectile  solution  data.  Table  II  lists  the 
parameters  of  the  nominal  firing  situation. 
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TABLE  II 

Example  Firing  Situation  Parameters 


Attacking  aircraft: 

Wing  angle  of  attack 

1.52  deg 

G-loac 

1.00 

Sight  depression  angle 

9.91  mils 

Sight  ceflection  angle 

0.00  mils 

Attituc'e  (stability  axes): 

Heading  (yaw) 

0  deg 

Pitch 

2.01  deg 

.  Roll 

0  deg 

Projectile: 

Muzzle  velocity 

3380  ft/sec 

Initial  velocity 

4331.38  ft/sec 

Initial  range 

2000  ft 

Final  velocity 

2915.97  ft/sec 

Final  range 

2630.12  ft 

Average  velocity 

3530.82  ft/sec 

Time  of  flight 

.7449  sec 

Calculation  of  Projectile  Miss  Distance 

If  the  aircraft  cannon  is  fired  when  the  attacking  aircraft  is  per¬ 
turbed  from  the  nominal  lead-pursuit  trajectory,  the  projectiles  will  miss 
the  target.  The  determination  of  terminal  miss  distance  is  rather  straight¬ 
forward  and  proceeds  in  the  following  manner. 

First,  the  aircraft  is  placed  on  the  nominal  trajectory  in  the  com¬ 
puter  algorithm  .and  then  it  is  perturbed  from  this  state.  The  disturbed 
state  is  generated  by  combining  the  nominal  aircraft  controlled  variables 
with  the  variances  of  the  perturbation  controlled  variables  determined  in 
Chapter  III  (Eqs  113-118);  however,  great  care  must  be  exercised  to  maintain 
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Then,  the  projectile  Is  fired  from  this  perturbed  state  and  gravi¬ 
tational  acceleration  and  aerodynamic  drag  fjrce  act  on  it  during  the 
time  of  flight.  The  target  terminal  position  is  determined  and  the  closest 
approach  of  the  projectile  is  calculated.  Finally,  the  total  miss  distance 
is  resolved  into  elevation  and  azimuth  miss. 

For  the  example  firing  situation,  the  miss  distance  computation 
yields 

elevation  miss  =  8.62  ft  (152) 

azimuth  miss  =  23.56  ft  (153) 

Equations  (152)  and  (153)  represent  the  completion  of  the  analytical 
portion  of  this  study.  The  entire  process  consists  of  disturbing  the 
attacking  aircraft  from  its  nominal  lead-pursuit  trajectory,  determining 
the  size  of  the  perturbations  in  the  controlled  variables  with  the  pilot 
actively  attempting  to  zero  these  errors,  and  finally,  propagating  these 
tracking  errors  into  terminal  miss  distance.  The  significance  and  sources 
of  the  projectile  miss  distance  will  be  discussed  in  the  next  chapter. 

This  method  of  determining  miss  distance  is  not  rigorous  and,  in 
effect,  uses  a  statistical  characteristic  of  random  variables,  the  variances 
of  the  controlled  variables,  as  deterministic  quantities.  There  is  no 
guarantee  that  the  correct  phasing  or  proper  magnitudes  of  the  perturbed 
variables  are  accurately  represented  by  this  procedure.  It  does,  however, 
give  a  qualitative  feel  for  terminal  miss  distance,  and  one  can  use  this 
data  to  decide  whether  terminal  error  is  significant  enough  to  warrant 
the  expenditure  of  time  and  money  to  correct  system  deficiencies.  More¬ 
over,  the  control led-variable  and  tracking  variable  intensities  do  provide 
an  accurate  measure  of  system  performance  and  can  be  used  to  evaluate  sys¬ 
tem  improvements. 
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To  correctly  determine  terminal  miss  distance  would  require  a  rather 
complex  simulation  of  the  example  firing  situation  and  this  task  could 
not  be  accomplished  within  the  limits  of  this  study. 
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V.  Conclusions  and  Recommendations 

Introduction 

In  this  copter,  the  results  of  the  analysis  will  be  evaluated  and 
then  compared  w  th  the  baseline  case  of  strafing  in  an  F-4C  (Ref  17: 
143-146).  Next,  recommendations  for  extensions  of  this  study  will  be 
given;  finally,  some  general  remarks  on  certain  mechanical  aspects  of 
the  analysis  thecry  will  be  made. 

Analysis  of  Results 

The  objective  of  this  study  was  to  develop  a  model  of  the  pilot- 
aircraft  system  which  related  the  pilot  tracking  performance  to  the  air¬ 
craft  dynamics.  This  was  accompl i shed  by  developing  two  mathematical 
models  which  demonstrated  the  propagation  of  a  disturbance  input  through 
the  pilot-aircraft  system  to  the  aircraft  controlled  variables. 

In  the  longitudinal  case,  the  analysis  showed  that  perturbation  pitch 
attitude,  the  longitudinal  tracking  variable,  had  a  root  mean  square  inten¬ 
sity  of  0.1308°  and  the  perturbation  angle  of  attack  had  an  rms  intensity 
of  0.1443°.  If  remnant  threshold  is  defined  to  be  the  rms  intensity  of 
pilot  controlled  output  in  the  absence  of  a  disturbance  input,  Eq  (74) 
can  be  used  to  calculate  a  threshold  value  of  .1118°.  It  then  becomes 
evident  that  remnant  accounts  for  approximately  90X  of  the  longitudinal 
tracking  error;  therefore,  it  is  of  questionable  value  to  alter  the  longi¬ 
tudinal  dynamics  to  achieve  only  small  increases  in  tracking  accuracy. 

In  the  lateral  case,  however,  Eqs  (115)  through  (118)  demonstrate 
that  remnant  accounts  for  only  20"  to  33%  of  the  rms  intensity  of  the 
controlled  variables  and,  thus,  the  improvement  of  lateral  dynamics  would 
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significantly  reduce  pilot  tracking  error  in  the  lateral  axis.  The  next 
question  that  arises  is  the  determination  of  the  proper  method  to  improve 
lateral  dynamics. 

The  most  obvious  method  of  correcting  ihe  deficiencies  of  the  lateral 
dynamics  is  to  perform  a  parameter  optimization  study  on  the  lateral  augmen¬ 
tation  systems  to  determine  the  system  gains  which  will  minimize  lateral 
tracking  error.  A  more  sophisticated  approach  would  be  to  investigate 
different  stability  augmentation  configurator  s  to  determine  the  feedback 
scheme  or  blend  of  schemes  which  affect  the  greatest  reduction  in  lateral 
tracking  error. 

In  Chapter  UI  it  was  determined  that  the  longitudinal  dynamics  con¬ 
strained  the  homogeneous  pilot  to  operate  in  a  limited  range  of  gains; 
this  resulted  from  the  effects  of  the  dynamics  of  the  fire  control  system, 
lo  illustrate  this  point,  the  fire  control  system  dynamics  were  neglected 
and  a  special  longitudinal  case  was  analyzed.  The  results  of  this  section 
indicate  that  the  lateral  pilot  model  parameters  dominate  the  homogeneous 
pilot  model  when  the  fire  control  system  is  disengaged.  However,  comparison 
of  Eq  (115)  through  (118)  with  Eq  (79)  through  (82)  shows  that  only  a  5 % 
improvement  in  lateral  controlled  variable  intensities  is  achieved  in  this 
case;  additionally,  Eq  (120)  indicates  that  the  rms  intensity  of  the  longi¬ 
tudinal  tracking  error  actually  increases  slightly  in  this  situation.  Thus, 
improvement  of  longitudinal  dynamics  will  not  materially  improve  lateral 
tracking  error. 

In  sjmmary,  the  F-4E  lateral  dynamics  introduce  significant  lateral 
tracking  errors  in  the  example  firing  situation  and  should  be  augmented 
to  improve  mission  performance.  Conversely,  longitudinal  tracking  error 
is  predominated  by  pilot  remnant  and,  thus,  improvements  in  aircraft 


GGC/EE/72-3 


dynamics  will  not  significantly  improve  longitudinal  tracking  performance. 
Baseline  Comparison 

Comparison  of  the  tracking  errors  in  the  aerial  gunnery  situation 
with  those  of  the  baseline  strafing  mission  will  give  an  indication  of 
the  capability  of  the  control  system  to  perform  different  missions  and 
will  il lustra t;e,.corutro-l  ■system  characteristic?  which  are  common  to  the 
two  missions.  The  longitudinal  tracking  error,  the  perturbation  pitch 
attitude,  is  larger  at  .26°  in  the  aerial  gunnery  case  than  in  the  strafing 
case  where  it  is  .17°  (Ref  16:133).  It  should  be  noted,  however,  that 
pilot  threshold  remnant  and  visual-motor  coupling  from  the  lateral  track¬ 
ing  task  dominate  these  longitudinal  tracking  error  values;  consequently, 
longitudinal  tracking  performance  cannot  be  significantly  improved  by 
.altering  the  pitch  tracking  dynamics. 

Similarly,  the  lateral  tracking  error  is  .41°  and  .16°  for  the 
aerial  gunnery  and  strafing  tasks,  respectively  (Ref  16:133).  In  this 
case,  the  lateral  tracking  error  in  aerial  gunnery  is  more  than  double 
that  in  strafing  and  is  due  primarily  to  the  differences  in  pilot  control 
pol icy. 

Although  the  tracking  errors  are  larger  in  both  instances  for  the 
aerial  gunnery  task,  the  same  dynamic  characteristic  manifests  itself 
in  both  missions:  lateral  tracking  error  predominates.  Hence,  if  lateral 
tracking  error  is  reduced  by  altering  the  lateral  tracking  dynamics,  longi¬ 
tudinal  tracking  error  will  also  be  reduced  as  a  result  of  the  cross  coupling 
effects.  A  design  goal  common  to  both  missions,  then,  is  the  improvement 
of  lateral  tracking  dynamics. 
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Recommendations  for  Further  Study 

The  purpose  of  this  study  was  to  demonstrate  the  effect  of  aircraft- 
dynamics  upon  pilot  tracking  capability.  As  evidenced  by  the  discussion 
in  Chapter  II,  however,  it  appears  that  computational  errors  in  the  fire 
control  system  may  contribute  significantly  tc  terminal  miss  distance. 

The  next  step  in  the  analysis  of  the  aerial  gunnery  mission  then  would 
•be  to  model  the  fire  control  system  and  demonstrate  its  effect  on  terminal 
miss  distance.  The  analyst  could  then  focus  his  attention  on  the  system-- 
flight  control  or  fire  control--which  had  the  greater  effect  upon  mission 
success. 

With  both  systems  completely  modeled,  the  sight  smoothing  could  be 
varied  to  determine  the  effect  upon  longitudinal  tracking  capability; 
also,  as  discussed  previously,  various  augmentation  schemes  could  be 
employed  to  reduce  lateral  tracking  error. 

Discussion  of  the  Analysis  Procedure  ' 

The  key  element  in  this  analysis  is  the  analytical  model  of  the 
human  pilot.  The  example  firing  situation  was  specifically  chosen  to 
meet  all  constraints  imposed  by  this  technique.  Nevertheless,  in  several 
instances  during  this  study,  the  pilot  modeling  theory  seemed  inadequate. 

In  both  lateral  and  longitudinal  axes,  system  dynamics  were  encountered 
which  were  much  more  complex  than  the  plant  dynamics  used  in  the  pilot 
modeling  studies.  Further,  in  the  longitudinal  axes,  the  pilot  was  con¬ 
strained  to  operate  in  a  narrow  frequency  range;  no  data  existed  to 
indicate  whether  the  pilot  could  really  operate  in  such  a  situation. 

Finally,  all  operating  frequencies  were  much  lower  than  those  in  the  studies' 
by  McRuer  (Ref  12) . 

Because  the  analysis  theory  used  in  this  study  has  so  many  potential 
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applications  and  because  the  pilot  modeling  techniques  are  such  a  key 
element  in  the  theory,  it  is  recommended  that  further  studies  using  more 
complex  dynamics  and  requiring  operations  at  lower  frequencies  be  con¬ 
ducted  tc.  expand  the  existing  body  of  knowledge  of  analytical  pilot 

modeling  theory 
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F-4E  Aircraft  Characteristics 

This  appendix  serves  as  a  compendium  of  information  of  the  charac¬ 
teristics  of  the  I-4E  aircraft  in  the  example  firing  situation.  The 
F-4E  can  be  mathenatically  described  by  its  dimensions,  its  control 
system  dynamics,  clid  the  coefficients  of  its  equations  of  motion  at  a 
particular  flight  condition.  The  flight  condition  of  interest  in  this 
study  is  the  nominal  lead-pursuit  trajectory  of  the  aerial  gunnery 
mission;  the  parameters  of  this  state  are 

h  =  15,000  ft  above  mean  sea  level  =  altitude 
M  =  .9  Mach  =  Mach  number 
Uq  =  951.63  ft/sec  =  velocity 
p  =  .001496  slugs/ft3  =  air  density  (Ref  6:497) 

The  dimensions  of  the  aircraft  at  this  flight  condition  are  (Refs 
3  and  11) 

W  =  46,159  pounds  =  aircraft  weight 

e  =  2  degrees  =  gun  depression  angle  below  the  fuselage  reference 
line 

Ix  =  28,203  slugs  ft2  =  moment  of  inertia  about  the  X  axis 
Iy  =  164,258  slugs  ft2  =  moment  of  inertia  about  the  Y  axis 
Iz  =  153,899  slug  ft2  =  moment  of  inertia  about  the-Z  axis 
Ixz  =  4427  slug  ft2  =  product  of  inertia 
c  =  16.04  ft  =  mean  aerodynamic  chord 
s  =  530  sq  ft  =  reference  wing  area 
b  =  38.41  ft  =  wing  span 
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The  dimensional  aerodynamic  coefficients  for  the  equations  of  motion 
written  in  stability  axes,  Eqs  (90)  and  (20),  are  listed  below  for  the 
example  firing  stfite.  The  longitudinal  axes  stability  coefficients  are 
(Ref  11:13) 


awj>  =  1*8  degrees  =  trim  wing  angle  of  attack 

C.G.  =  31.9%  =  center  of  gravity  expressed  as  a  percentage  of 
•mean  aerodynamic  chord 

Za  =  -  1.0326 

Zg  =  -.005291 

Z5  =  -.09512 

^  =  -.3439 

Mo  =  -10.443 

Mg  =  -.7381 

M6  =  -37.08 

The  lateral -directional  axes  stability  coefficients  are  (Ref  3:24) 


Lp  =  -  1.972 
4  =  .5334 
4  =  -  21.00 
L6r  =4.102 
L6  =  13.57 
Np  =  -  .0314 
Nr  =  -  .4674 
N0  =  10.73 


Yp  =  .0002 
Yr  =  -  .9982 
Y6  =  -  .1886 
Y6r  =  .0250 
Y6  =  -  .0040 
Y^  =  .0338 
N$r  =  -  5.039 
N6  =  .3857 


Ic  should  be  noted  that  the  subscript  6  indicates  the  combined  effective 
displacement  of  spoilers  and  ailerons  for  the  F-4E;  in  this  study  it  will 
be  referred  to  simply  as  aileron  deflection. 
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Fig.  23.  F-4E  Elevon  Control  System  (Adapted  from  Ref  17:155) 
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The  aileron  control  system  of  the  F-4E  is  described  in  Fig.  25 
(Ref  17:157). 


POWER 

CYLINDER 


Fig.  25.  F-4E  Aileron  Control  System  (Adapted  from  Ref  17:157) 

The  control  law  for  the  system  is 

6  _  10 5 PA  _  2.65p 

a~s  +  10~s  +  10 

This  can  be  expressed  as 


<5  =  ACsltS^.Cs)  -  B(s)d] 
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